


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations l. Thesis and Dissertation Collection, all items 


1956 


Helium gas turbines for nuclear power. 


Phillips, Kenneth Elwin 


Monterey, California: U.S. Naval Postgraduate School 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
۳ W D U DLEY research materials and institutional publications created by the NPS community. 
FW i Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


1 ۱ KNOX appointed — and published — scholarly author. 
| In D 
| LIBRARY Dudley Knox Library / Naval Postgraduate School 


411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


HELIUM GAS TURBINES FOR NUCLEAR POWER 





Kenneth Elwin Phillips 


Library 
U. S. Naval Postgraduate School 
Monterey. California 























HELIUM GAS TURBINES 
FOR 


NUCLEAR POWER 


Kenneth E. Phillips 





NUCLEAR POWER 


by 
Kenneth Elwin Phillips 


Lieutenant, United States Navy 





Submitted in partial fulfillment 
of the requirements 
for the degree of 
MASTER OF SCIENCE 
IN 
MECHANICAL ENGINEERING 


United States Naval Postgraduate School 
Monterey, California 


6 و ۶ ۲ 





PREFACE 


In recent years great progress has been made in the gas turbine 
field. The improvements made in high temperature materials, compressor 
and turbine design, and compact heat exchangers have made the gas tur- 
bine power plant a serious competitor of steam plants and diesel engines 
for many applications. There are, however, many problems still exist- 
ing that tend to make conventional gas turbines unattractive. 

The use of residual fuels in the open cycle results in ash deposits 
and corrosion in the turbines which seriously limits the operating life 
of the units. Furthermore, open cycle plants are limited in power rat- 
ing since the components become too large for outputs greater than about 
20,000 HP. The closed cycle gas turbine plant eliminates these problems 
since only clean air circulates through the cycle and high pressures 
permit components of a practical size for outputs much larger than 
20,000 HP. For a ship propulsion unit where the power requirements 
vary over a wide range, the closed cycle has the additional advantage 
of good part load efficiency. A major factor in the slow progress of 
the closed cycle plant has been the inability to manufacture a reliable 
air heater cheaply enough to permit the gas turbine plant to compete 
with conventional steam plants. 

The future of the gas turbines may depend on how well they can be 
adapted for use with nuclear power. At the present time the only ship 
propulsion plant employing a nuclear reactor is essentially a steam 
plant. The size of this plant is large since a two loop system is 


required in order that energy may be conveyed from the reactor to the 
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secondary working fluid. A simpler and smaller plant would be one in 
which there is only one fluid serving both as the working fluid and 
reactor coolant. Such a plant is a closed cycle gas turbine plant 
using helium at high pressure. Helium has good coolant characteristics 
including a negligible cross section for neutron absorption. If pure, 
it is not subject to induced activity. By proper design of the cooling 
channels and filtering system, fission products can be prevented from 
being carried through that portion of the cycle which is external to 
the reactor. This would require only the reactor to be shielded. It 
has been presumed up to this point that helium is a satisfactory gas 
turbine fluid. Although helium turbines have never been built, the 
following general statements can be made by comparing the properties of 
helium with air: 

(1) Better heat transfer characteristics of helium results in 

smaller heat exchangers. 

(2) Lower molecular weight of helium requires more stages in the 

turbo-machines for the same peripheral speed. 

(3) Higher velocity of sound in helium permits higher peripheral 

speeds before Mach effect occurs. 

It would appear from these statements that the advantage of smaller 
heat exchangers for a helium cycle may be more than offset by requiring 
more expensive turbines and compressors as well as a costlier and less 
available gas. Design criteria for helium compressors are T on 
stress limitations, whereas for air compressors the Mach effect becomes 
the limiting factor. Thus, the higher peripheral speed permitted for 


the helium compressors will partially counteract the effect of lower 
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molecular weight on the number of stages required. Furthermore, improv- 
ed efficiency and fewer stages of compression may result by employing 
a higher effectiveness regenerator than would be practical with air. 

It is the writer's objective to investigate a 50,000 HP closed- 
cycle helium gas turbine power plant for ship propulsion and determine 
the approximate sizes of the major components. Unfortunately, the 
scope of this report does not cover many of the interesting problems 
which would have to be investigated and solved before this plant could 
become a reality. Some of these are discussed briefly in Chapter IX. 

It has been necessary to include a list of symbols at the begin- 
ning of each chapter since the same symbols may have different meanings 
in different chapters. 

The writer wishes to express nis deep appreciation for the encour- 
agement "neon و‎ given him by Professor M. H. Vavra and Professor 
H. M. Wright in the preparation of this report. The writer is also 
grateful for the fine clerical support given him by his typist, Mrs. 


Marílyn Heard. 
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CHAPTER I 


THERMODYNAMIC PROPERTIES OF HELIUM 


Symbols used in this chapter: 


S ol Description 


B. virial coefficient 
E Specific heat at constant pressure 
Cy specific heat at constant volume 
f reheat factor for OO stages 
h enthalpy 
k thermal conductivity 
P pressure 
R gas constant 
T absolute temperature 
t temperature 
V specífic volume 
Y ratio of Cp 

c 

V 
A increment 
7) efficiency 
H absolute viscosity 
Y corresponds to Y -1 

Y 

Subscript Definition 
0 reference quantity 
C compressor or compression process 


Units 
£t? /1b 
BTU/lb - °F 


BTU/lb - °F 


a 


BTU/ft^ - hr - °F/ft 


psia 


£t^/1b 


lb/hr - ft 





Subscript Definition 


i initial 

E final 

S isentropic process 

t turbine or expansion 
General 


In studying a cycle that uses fluids such as steam or air, or one 
that behaves like an ideal gas, there are tables, charts and equations 
which can be used to obtain accurate thermodynamic information. This 
is not the case for helium. Several equations of state have been 
proposed based on theory and empirical data obtained at low tempera- 
tures. Which of tnese equations, if any, is applicable at higher 
temperatures will be known only after more extensive empirical data 
are available. Until then, extrapolation from the lower temperatures 
will have to suffice. The extrapolation methods used in this investi- 
gation are based on the research done by Akin (1) His tables for 
enthalpy, entropy, and specific volume extend to 600ºF whereas the 
highest temperature in the cycle being investigated is 1350ºF. His 
extrapolation formulas for absolute viscosity and thermal conductivity 


extend to temperatures beyond 1350°F. 


Absolute Viscosity, 4 
(2) 


Akin found Keesom's equation for 4 accurate within 1% for 


temperatures between -450 and 2000°F. This equation, plotted in 
Figure 1-1, is: 


u = 3,315 er x rome [b/hr-ff (1-1) 


Thermal Conductivity, k 


The following equation is based on kinetic theory and test results 


at temperatures below 212ºF. It is plotted in Figure 1-1. 


K= هر‎ 6۷ 6 SNUG le = a 
or k= 18/66 ۸ Duo E (1-2) 
Specific Volume, v 
Akin's results on specific volume are based on the Beattie-Bridge- 
man equation of state. Comparisons with the co-volume equation of 
state show good agreement in the range of pressures and temperatures 
encountered in this report. Because of its simpler form, the latter 


equation has been used throughout. 


P= m psia (1239 


d 3 
where B.= 0.306 (= F - 1845 (4% 0.822 (+) (1-3a) 


NO 


and R= 2.6829 ft Ibe 
En 

Equation (1-3a) is plotted in Figure 1-2. The specific volume for 
temperatures from 0ºF to 1400ºF when calculated from equation (1-3) 
is 0.047 to 0.040 £t? /1b greater than that obtained by assuming ideal 
gas behavior at the same temperature. 

Enthalpy, h 

For enthalpies below 600?F, Akin's tables have been used. The 


extrapolation formula for a constant pressure process can be expressed 


— 





le Cs (t-t) + ۲ ) رط‎ P) SER 
where h, = h at - 440 F 
and Co = JE BTU/Ib- F 


At temperatures above 400°F it is apparent from the tables that for a 


constant pressure process: 


h= h, + c (tt) BTU (1-5) 


where t > 400 FF 
ad Cp = h24 BTU b-'F 


Since extrapolation is necessary above 600ºF, (1-5) becomes: 


he E 7 (t - 600) BTW lb (1-53) 


Figure 1-3 contains plots of h vs. P for temperatures from 100ºF to 


600°F. 


Isentropic Work of Expansion and Compression 
For an ideal gas the work of compression for a reversible-adiaba- 


tic process is: 


2 ۰ F BTU/lb (1-6) 








Similarly, the work of expansion is: 
Y-! 


apa Ties 2 ۱ BTU/ Ib (1-7) 


An analysis of Akin's tables reveals that these equations are not 
applicable for helium if both > and Y are assumed to be constant. 


Since p. is constant (1.24 BTU/1b °F) equations (1-6) and (1-7) were 





E using various combinations of Ti, Pi, and MP 


Y 
solved for = 





Y - 
The results indicated that 


we is a function of Pi only. Since this 


En 


function differs for the compression and expansion processes, X 


has been replaced by Ye and M respectively. 


- ۸ 
Ah. = n E S z BTU/Ib (1-6a) 


e 55 - (2) BTU/|b (1-7a) 


Figure 1-4 contains plots of bg and ae VS. Poe 


4 


Reheat Factor, l+f ao 
The values of ae and y, are also applicable to expressions for 
evaluating reheat factors. Assuming ideal gas behavior and replacing 


Y -i 
E — by Mee or Y, the reheat factor, l+fo9 , becomes: 


W 
4 
25 « 
A al es SER for compression 


(a) 
(2 — | (1-8) 


If. Es N 


in 





h A A for expansion 
t 


5 
| — (2) E (1-9) 


A more detailed analysis of reheat factors is given in Chapter IV on 
Compressors", 
By means of the aforementioned formulas and/or the plots, all the 


thermodynamic properties that are required in the following cycle 


analysis can be obtained. 





OT X (33/4, - ,33 - 1u/nId) 'x 'Kar^4T39npuo) ] 0 
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CHAPTER 11 


PRELIMINARY CYCLE ANALYSIS 


Symbols used in this chapter: 


Symbol 


STE ات‎ 


Description 


virial coefficient 


specific heat at, constant pressure 


total isentropic enthalpy change for 
the compressor or turbine 


enthalpy 

thermal conductivity 
mass rate of flow 
pressure 


gas constant 


absolute temperature 
temperature 
specific volume 


ratio of "p 


C 
V 


increment 
efficiency 


absolute viscosity 


IT 


corresponds to = 





specific heat at constant volume 
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Units 
FE /ID 
BTU/lb - ?F 


DTU/lb - °F 


BTU/1b 

BTU/1b 

BTU/£t? - hr - *F/ft 
1b/sec 

psia 


ft - Ib, 


lb - °F 
m 


°R 
°F 
tt^ /1b 


m de m es 


lb/hr - ft 





Subscript Definition 


0 reference quantity 
e compressor or compression process 
i initial 
E final 
R regenerator 
S isentropic process 
t turbine or expansion 
Th A thermal 
General 


Since it is the objective of this report to analyze a 50,000 HP 
closed-cycle a turbine plant, it is of prime importance to select a 
set of conditions which will result in a plant of high thermal effi- 
ciency and reasonable size which is still practical and economical to 
build, operate and maintain. To make such a selection one must know 
the following: 

(a) Components of the cycle 

(b) Reasonable estimate of the machine efficiencies 

(c) Anticipated losses associated with the various components and 

interconnecting piping. 

(d) Regenerator effectiveness 

(e) Limiting temperatures, both upper and lower 

(£) Highest pressure permissible in the cycle. 

Components of the Cycle 


This is not a study of various types and arrangements of cycles, 


hence a cycle of minimum complexity was chosen. The basic cycle is 


13 


similar to the conventional plant built and operated by Escher-Wyss (3) 


and consists of a single heat source, regenerator, precooler, compress- 
ors with intercooling, and two turbines. For the heat source, the 
conventional air heater is replaced by a homeogeneous reactor, the core 
consisting of a graphite and highly enriched uranium mixture. Since 
auxiliary equipment, including control apparatus, has little effect on 
the thermodynamics of the cycle it is not included here. A diagram- 
mee sketch of the basic cycle is shown in Figure 2-1. 

Machine Efficiencies 

Assigning efficiencies to the two compressors and two turbines 
can at best be only a first approximation. Turbo-machines are built 
today with overall efficiencies over 80% when used in conventional 
cycles, However, for a given cycle such as the one being inve ie 
practical design consideratíons may preclude these hígh efficiencies 
and a conservative estimate is in order. No distinction is made in 
the preliminary analysis between internal and overall machine effi- 
ciencies and no separate mechanical efficiency is specified. The 
efficiencies for Compressor No. 2 and Turbine No. 2 are lower than 
for the No. l units in the event that radial flow turbo-machines are 


practical. 


MACHINE EFFICIENCY (Assumed) 
Turbine No. 1 0.85 
Turbine No. 2 0.80 
Compressor No. 1 0.84 
Compressor No. 2 0.80 
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Pressure Drops 


The losses associated with each component are expressed here as a 
ratio of the static pressure loss across the component to the inlet static 
pressure. Again, these are only estimates based on rough calculations 
and experience and include not only the drop across the component but 


also associated inlet and exit piping effects. 


COMPONENT a (Assumed ) 
i 

Reactor 0.010 
Regenerator 

High Pressure Side 0:015 

Low Pressure Side 0.010 
Precooler 0.0075 
Intercooler 0.005 


Regenerator Effectiveness 


The regenerator effectiveness is defined as: 


no E ول‎ h4 )2-1( 
R hg -h4 

A high effectiveness is usually desirable from the standpoint of high 
thermal ef.iciency. It is not practically feasible to have De greater 
than 0.9 because of the large amount of surface required. In the 
preliminary analysis N R is taken as 0.7, 0.8 and 0.9 to indicate its 
effect on efficiency. 

Limiting Temperatures of Working Fluid 

The lower temperature is dependent primarily on the cooling medium, 


in this case sea water. For a given plant a decrease in the lower 





temperature will improve plant efficiency. However, since a ship may 
operate in tropic waters as well as arctic waters, a conservative 
lower limit of 100ºF is chosen. The upper temperature limit depends 
on materials and is gradually increasing as materials and cooling 
techniques are improved. Temperatures of 1450 and 1500ºF are common 
in present-day gas turbine cycles, but because of the high pressure 
encountered in this closed cycle an upper temperature of 1350°F is 
assumed. 

Maximum Pressure 

The maximum pressure in the cycle has a slight effect on effi- 
ciency due to the variation of M. and pat with pressure. More 
important, however, is the effect of pressure on the size of the plant. 
Higher pressures mean smaller specific volumes, hence smaller flow 
areas. Based on steam plant design a maximum cycle pressure of 1000 
psia seems reasonable for a plant of this size. 

Cycle Analysis 

Now, with the foregoing assumptions, together with the thermody- 
namic information from Chapter I, the cycle of Figure 2-1 has been 
analyzed by standard methods to obtain thermodynamic properties at 
each point for pressure ratios between 1.25 and 5.0. In addition, 
the following performance indices are plotted in Figures 2-2, 2-3 and 
2-4 respectively: 

(a) Thermal efficiency based on work output at turbine shaft 

(b) Helium rate 


(c) Work ratio 
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A list of equations used in the analysis 


V = 2,6829 T 7 D. 


AH = Er 


h, = heor + 124(t-600) BTU Ib 


follows: 


Hb 


o E 
AH,— 12% Ti E 1 BTU/b 


ln 
AH = L24 T ۱ - (E) oro 


h = fuhe 
th he-hs 





Helium Rate = 2545 == 
Helium Rate = : Ib / hp hr 


(1-3) 


(2-2) 


(2-3) 


(1-5a) 


(1-6a) 


(1-72) 


(2-4) 


(2-5) 





Work Ratio = hi-hs (2-6) 


E ie he 


n — hz-ha (290) 
fs n 
8 + 
Conclusion 


In studying Figure 2-2 we observe the following points of maximum 


efficiency: 
P 
AR VM er 
0.7 3.00 0.268 
0.8 2.65 0.292 
0.9 2.25 0.329 


If this cycle is to compete with other power plants it must have a 
higher thermal efficiency and/or be smaller in size. A high effective- 
ness regenerator provides the higher efficiency at the expense of 
increased size of plant. At the lower effectiveness, not only is the 
thermal efficiency less, but the larger pressure ratio will result in 
several more stages for the turbo-machines. Placing emphasis on the 
higher efficiency and fewer stages, the remainder of this report consid- 
ers only the cycle having "up, 2 wand fi = 0.9. Table 2-2 


R 


lists the thermodynamic properties at each point in this cycle. 
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TABLE 2-1 


PLANT PERFORMANCE AS A FUNCTION OF PRESSURE RATIO 


Pressure Plant Thermal Efficiency, H th Helium Work 
Ratio Rate Ratio 
"Alp, "& -0.9 "]n = 0.8 "R - 0.7 (1b/hp-hr) 

1:25 0.184 0:122 0.091 S10 0.378 
1.50 0.279 0.210 0.169 24.4 0.413 
[75 0591222 * 0.253 0212 Ln 0.410 
2.00 0.323 0.274 0257 14.8 0.400 
2525 0.329 05285 0.252 13:2 0.388 

2.50 0.327 0.290 0.261 1202 02375 
2957 5 0.323 0.291 0.265 11.6 07562 
3.00 0.318 0.291 0.268 LIO 0,381 
4.00 0.290 0:275 0.261 10.4 0.308 
5.00 0.261 0.252 0.244 10.6 0.270 
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TABLE 2-2 


THE RMODYNAMIC PROPERTIES FOR CYCLE HAVING LA 2 2.25 AND Pn 21099 


Location P (psia) h(BTU/1b) t( °F) vere IB) 
(Fig.2-1) 
1 444.4 701.5 100 3.428 
2 666.7 648.8 216.9 2.768 
3 663.4 703.7 100 2:309 
& 1000 861.9 224.9 1.883 
3 985 1667.8 874.8 3.677 
6 27122 2256.9 1350 22019 
7 631,2 1951.0 1106.3 6.692 
8 452.2 1757.3 ge 8.413 
9 447.7 951.4 »01.> 4.607 
m x 182.4 lb /sec (For Net output = 50,000 HP) 
Nth = 0.329 
Helium 
Rate = 159 lb/hp-hr 
Work 
Ratio = 0.388 
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CHAPTER III 


HEAT EXCHANGERS 


Symbols used in this chapter: 


Symbol 


A 
Q 
E 


3. 


Nu 


Pr 


Re 


Description 


heat transfer area per tube 

free flow area per tube 

frontal area per tube 

fin longitudinal length 

fin thickness 

specific heat at constant pressure 
heat exchanger shell inside diameter 
hydraulic diameter 

tube diameter 

convection heat transfer coefficient 
thermal conductivity 

heat exchanger tube length 

fin radial length 

(defined by equation 3-10a, page 33) 
mass rate of flow 

total number of tubes 

Nusselt number 

pressure 

Prandti number 

Reynolds number 


tube center to center spacing 
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Units 

ft 

ft 

ft 

in. 

in. 

BTU/lb - °F 

ft 

ft 

in. 

BTU/£t* - hr - ^F 
BTU/£t* - hr - °F/ft 
ft 

in. 

lb/sec 

psia 


in. 





PS Ss 


Subscript 


C 


Description 


temperature 
overall heat transfer coefficient 
velocity i 
specific volume 
ratio S RE y 

(ie), 


wall thickness or boundary length 
thickness 


increment 

friction factor 

fin efficiency 

fin temperature effectiveness 
absolute viscosity 

density 


dimensionless overall heat transfer 
coefficient 


dimensionless convection heat 
transfer coefficient 
A 


۹ 


ratio of heat transfer areas, 


ratio of fluid flow areas, 


SIS 
xc 


Definition 
cold fluid 
exit 
fin 


hot fluid 
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Units 


°F 

2 8 
BTU/ft” - hr - °F 
ft/sec 


FE IND 


"B ou «mo 


in. 
lb/hr - ft 


IUE. 


“m «ue VP oom 


æ Ge az 





| 


nn gg ہہ‎ 





Subscript Definition 


i inlet, inside, or initial 
0 outside 
t tube 
W wall 
General 


Having selected a pressure ratio of 2.25 from the preliminary 
analysis, the components will now be more closely investigated to 
determine if they can be of practical design. Since one of the more 
important requirements is the size of the plant, a design study of the 
regenerator, intercooler and precooler is necessary. Before going into 
this study, some review of heat transfer theory as used herein may be 
advantageous. 

Heat Transfer Theory 

In the usual type of heat exchanger heat transfer occurs between 
two fluids of different temperatures separated by a wall of finite 
thickness. The quantity of heat, dQ, transmitted through an element of 


area, dA, can be expressed as: 


SR = U(t,- td w 


where U is the overall heat transfer coefficient associated with are: 
dA. 


lf the wall surface is a plane, i.e., no extended surface: 


2 7 O + ot (3-2) 
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C ال‎ N 
U, => h, + A ۰ "E Ac : (3-2a) 
An END É 


Kays has modified (3-2a) to include extended surfaces on both 


sides of the wall: 








== + K (329 
U A Aw A 
ha, í Mw à A, Noche 


where den A, = effective surface area on hot side 


and o. A. z effective surface area on cold side. 
For thin walled surfaces the term En is negligible: 
wW 


pk 


An 


l ۱ | | 
— == ct ل‎ (3-3a) 
Un Noy, N, Ag m h 
. Ay c'ic 


Combining the techniques of Kays and I it is convenient 


to change (3-3a) to a dimensionless form: 


E |] $‏ 
)3-4( کک + — حه سل 

eee -4a 
where 9 == PGV (3-42) 
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A h (also known as 





PV Stanton Number)  (3-4b) 

W = m (3-4c) 
P), 

D "T (3-4d) 

Li E (3-4e) 


Considering the hot fluid only, the heat loss, dQ, for a tempera- 


ture change dt can be expressed as: 


dQ, =- ce dt, = —-(eQve) dt, (355) 


Combining (3-5) with (3-1): 


Alp aes‏ - و94 
Pr Ay Cta‏ 





If we assume that the configuration and dimensions of the heat 
exchanger elements do not change along the direction of the flow, then 


dA, = ^n dx, where x is a lengthwise coordinate, L is length of the 
L 


surface and AL is the total hot side area. Integrating equation (3-6): 





L te 
Ph Ah ۲ dt (3-7) 
7 











or 





T 
RR ge 
à, fa d(2) | DA 8e (3-7a) 


Ph will change as the fluid flow proceeds along the surface and 
its function of x must be known before the integral can be evaluated. 
Likewise, tort, changes and its function of t must be known. 


Usually these functions are not known and it is common practice to 


simplify (3-7) as follows; 


An 1 = (ti- te), (3-7b) 
AT 


L 
where $, = / 4a4 





| 
dt 
toe 


Bowman, Mueller and Nag1e 9) have determined relations for At 


EP: 
and At = ua 
€ 


for many types of heat exchangers. For the longitudinal counter-flow 


exchanger At is the familiar log mean temperature difference: 


A (tna = teal = (tne e) (3-8) 
In ih ate 
(the- te 


Close approximation is obtained if $ is evaluated by equation 


(3-4) using average values of W and VW n 
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Up to this point the type of extended surface has not been speci- 
fied. Once this is specified there remains the problem of evaluating 
ns 

DE À has been defined as the effective heat transfer area. For a 
finned surface this effective area can be expressed as the sum of the 
wall area A, (not including fin base), and an effective fin area. 
From the definition of fin efficiency, De this effective fin area is 


p A; where Ar is the total surface area of the fin. 


A = Aw + ۷۶ ۸ 
re =A- Pi 


e ANS ۳ (i7 f) (3-9) 


For several types of fins the solutions for N are known. 
ran) 
As an example, for a bar fin of constant thickness shown in 
Figure 3-1: 
A 
toy 
e 


Figure 3-1 


n - fons A (3-10) 
m 


where m = | 2h(L+c) 
KELE 
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tE Le C e 2 (3-10a) 
Ke C 


It is noted that Gardner's solutions for h ç are based on several 
assumptions. Of these, the following may result in appreciable error: 
(1) The heat transfer coefficient, h, is uniform over the 
entire fin surface, 
(2) The temperature of the base of the fin is uniform. 
(3) Temperature gradients normal to the fin surface are 
neglected. 
(4) Heat transfer through the outermost end is negligible. 
The only type of heat exchanger involved in this investigation is 


tubular with external fins. By introducing the hydraulic diameter, Da: 


D = 4 Ch (3-11) 
WIL 
equation (3-7b) becomes: 


m | (ti zug (3-12) 


D, 48 AT 
The above equation will be one of the principal equations employed 
in this investigation for design of the heat exchangers. 
Pressure Drops 
Fluid pressure drop in the heat exchanger must be considered in 
addition to heat transfer aspects. It has been found both experimen- 


tally and analytically that the two are interdependent. Improvement in 
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heat transfer characteristics results in increased pressure drop, so a 
compromise has to be made in most designs. The relation used by 
Be MEAGRE S) and others for longitudinal flow inside of tubes 


is: 





AP = y E Y (3-13) 


where 5 is the coefficient of resistance. It is identical to 4f 
in McAdams. 
(8). 


For turbulent flow through smooth tubes 


0.2 


C= 0.184 Re. (3-14) 


Sooo < Re < 500000 


Similarly, for determining pressure drop outside of the tubes (on 
finned side), equation (3-13) and (3-14) are applicable if the diameter 


is replaced by an equivalent diameter: 


Do = 4 x free Flow area بح‎ 4h (3-11a) 
eq wetted perimeter Ay, 


h 


perimeter of the exchanger shell. The error in assuming D, pm" is 


The difference between D, and MEC is that Deg includes the wetted 


small for most heat exchangers. 
In addition to the pressure drops due to pipe friction there are 
Static pressure changes due to acceleration or deceleration of the flow 


at the heat exchanger ends and changes in specific volume. Since these 
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effects depend on pipe sizes, types of headers, etc., they are left for 


further discussion in Chapter VII. 


Correlation of Heat Transfer and Pipe Friction 


Equation (3-4b) can be expressed in terms of other dimensionless 


quantities as follows: 


۳ Nu (3-15) 





where Nu — ha (Nusselt No.) (3-15a) 
Pr == "m (Prandtl No.) (3-15b) 
dV 


(Reynolds No.) (3-15c) 


Re 


I 
| 


< 
< 


(8), 


For turbulent flow inside of tubes according to McAdams 


8 
Nu = 0.023 Re" Pr! (3-16) 


Equation (3-16) can be used for turbulent flow outside of tubes 
providing that the hydraulic diameter is used in lieu of d. 


From equations (3-14), (3-15) and (3-16): 


(3-17) 


To summarize, the principal working equations for tubular heat 





exchangers having longitudinal fins are: 





(3:995 


AP = $ de ME (3-13) 


Type of Heat Exchanger 


The feasibility of this cycle depends to a large extent on the size 
of the three heat exchangers; namely, regenerator, precooler and inter- 
cooler. By far the largest of these will be the regenerator. For an 
ordinary shell and circular tube regenerator analysis reveals that to 
satisfy the pressure drop and size limitations extremely small tubes 
(less than 0.2 in. diameter) placed very close together (less than 0.3 
in. center to center) are required. 

Construction difficulties, if the above dimensions were specified, 
would be of sufficient severity to render the design impractical.  How- 
ever, the addition of longitudinal fins to the exterior walls of the 
tubes in this type of heat exchanger will cause a salutory effect. 
Analysis reveals that tube size and spacing may be increased, and the 
total number of tubes decreased. Although this may not be as good a 
design as could be obtained by employment of more compact surfaces such 
(4) 


as were investigated by Kays, London and Johnson the reduced complex- 


ity in design and analysis warrants its use here. 


A similar design is investigated for the precooler and intercooler, 














the difference being in dimensions. In all cases, the hotter fluid is 
outside the tubes (on the finned side). Figures 3-3 and 3-4 show the 
finned tubes and tube arrangement respectively. Although the analysis 
assumes that the fins extend the entire length, L, of the tube, boundary 
layer buildup between the fins will require the use of several shorter 
fins staggered along the tube. This has been discussed by Glover m 


with the conclusion that the minimum fin spacing should be not less 


than one boundary layer thickness, 


7 





Turbulent Boundary 
Layer 


Figure 3-2 


Von ua) equation for turbulent boundary layer thickness, 


Ô» is: 
^s [vw 5 
S e uus A in, (3-18) 


max 
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Heat Exchanger Tube 


Figure 3-3 











Tube Arrangement 


Figure 3-4 
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For moderate Re: 


Vay = a (3-19) 
V 


Vmax > E Vow (3-194) 


In all three heat exchangers the tubes will have six longitudinal 


BTU 
re  {) 


: tt? - hr - °F/ft 


must have good metal-to-metal contact vith the tube. The assembly 


, 


brass fins (k,2 60 For good heat transfer the fins 
shown in Figure 3-3 consists of two concentric tubes, the outer one hav- 
ing the six external fins. These outer tubes will be quite short (six 
to ten inches) and could be made in one machine shop operation. These 
would then have to be brazed to the inner tube in a manner to insure 
intimate contact. 

In this analysis it has been found convenient to employ areas 
associated with one tube. For the equi-angular arrangement that has 
been selected, the frontal area per tube, Ol, is equal to twice the 


area of one of the triangles formed by the centers of three adjacent 


tubes. The formulas for this and other areas are: 





2 
کک‎ « e(24 -c)] 5 A (3-20) 


۱ 


۱2 4 E o 1 ft (3-22) 


= 
O سس‎ o m-—— 


a Tt 





rp» [> rp» 
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s? Sin 60° ne 
j4 | 


Q 9 bs E: + 6e | Ft 
h t 4 144 





Other equations which are used in the analysis are: 


Ve — Oy ve me (3-26) 
۷ Q. Yo, mn 
2 2 
D — 40eN ft (3-27) 
TT 
MY _. Meke (3-28) 
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Regenerator 


NS AN AN 
۸ an En 
T 2 \ 2 N 2 NEL 
/ / 
X N / 
/ 
\ , N / 
S XN 7 
SA NY 
2 8 
Figure 3-5 Schematic Representation of the Regenerator 
TABLE 3-1 


PROPERTIES OF HELIUM ENTERING AND LEAVING REGENERATOR 





E 3-5) t (^F) P (psia) h (BTU/1b) v (ft /1b) 
4 224.9 1000 361.9 1.883 
| 5 874.8 985 1667.8 7 
8 951.3 45212 1757.3 245 
3 ors R 951.4 4.606 
ES 


76.6 °F 


4 
it 
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Since this analysis is based on average values of such quantities 
as v, C? k, M » Re, Pr, Y , and $ , they are evaluated at the 
arithmetic mean temperatures. To distinguish the hotter fluid outside 
the tubes from the cooler fluid inside the tubes the terms "hot side" 


and "cold side" will be used henceforth. 


TABLE 3-2 


AVERAGE PROPERTIES OF FLUIDS IN THE REGENERATOR 


Quantity Hot Side Cold Side Units 
t 626.4 249.8 °F 
V 6.510 2.780 tt^ /1b 
k 0.138 0.132 BTU/£t” - hr - °F/ft 
JA 020759 020725 lb/hr - ft 
Pr 0.681 0.680 ---- 
E 1.24 1.24 BTU/lb - °F 
Precooler 


The basic difference between the regenerator and the precooler is 
that the fluid inside the tubes for the latter is sea water rather than 
helium. The heat transmitted to the water is lost, but the lower tem- 
perature for compression results in less compressor work required, and 
because of the higher temperature differential in the regenerator, in a 
higher overall cycle efficiency. The precooler, as in the case of the 
regenerator, consists of externally finned circular tubes with single 
pass counterflow. The tubes and tube arrangement are also the same as 


in the regenerator with a slight modification of dimensions. 
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Figure 3-6 Schematic Representation of the Precooler 


TABLE 3-3 


PROPERTIES OF FLUIDS ENTERING AND LEAVING PRECOOLER 


Location 


1 
9 
H,0 in 


H,O out 


Quantity 


E 


V 


Pr 


t (°F) P (psia h (BTU/1b) v (LE /1b) 
100 444 ,4 701.5 3.4208 
301.5 447.7 951.4 4.607 

60 m See 01604 
90 = ---- .01610 
à. = 182.4 1b/sec É 


At = 103 °F 


TABLE 3-4 


AVERAGE PROPERTIES OF FLUIDS IN THE PRECOOLER 


H,O 


Helium Units 
200.8 75 °F 
4.018 0.01607 et” /1b 
0.0558 2.22 ۱/۲ 6 
1.24 1.0 BTU/1b - °F 
0.10 0.346 BTU/ft^ - hr - ?F/ft 
0.692 6.42 en 
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Intercooler 





Figure 3-7 Schematic Representation of the Intercooler 


PROPERTIES OF FLUIDS ENTERING AND LEAVING INTERCOOLER 


Location t (CF) P (psia) h (BTU/1b) V (£t” /1b) 
2 216.9 666.7 848.8 2.768 
3 100 663.4 703.7 2.309 
H,O in 60 === ---- 0.01604 
11,0 out 90 ---- ---- 0.01610 
m = 182.4 lb/sec 
At = 75.4 °F 
TABLE 3-6 


AVERAGE PROPERTIES OF FLUIDS IN THE INTERCOOLER 
H_O 


Quantity Helium e Units 
t 15075 75 ?F 
V 2.538 0.01607 tt^ /1b 
A 050532 2,22 lb/hr - ft 
e 1.24 1.0 BTU/lb - ”F 
k 0.096 0.346 BTU/ft^ - hr - ^F/ft 
Pr 0.687 6.42 suec 
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Selection of Tube Dimensions 

The dimensions that appear in Table 3-7 are the result of many 
trial and error calculations. There are, no doubt, other combinations 
which will produce equally good results. It is intended that the ones 
used here represent a reasonable design for this cycle. 

Calculations 

Values in Tables 3-8 and 3-9 were obtained from substitution of 
known and assumed quantities into equations outlined earlier in this 
chapter. It can be seen that there occur several relations and vari- 
ables. However, there is only one independent variable and it is 
convenient to choose V as this variable. Figures 3-8, 3-9, and 3-10 


h 


show the variation of AP, L, D and total volume with V for the three 


units. 


۱ 


The final selection of V "depends on how large a pressure drop 


h 


and/or size can be tolerated. The choices of V, and the corresponding 


h 
heat exchanger dimensions listed in Table 3-11 are based on the values 
assumed for AP/P, in Chapter II. A knowledge of end effects was 
necessary to make the selection since in some cases the AP due to 


the end effects is larger than AP due to pipe friction. This is 


covered in more detail in Chapter VII. 
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ls 


ch T ea 


۲ o 
TUBE DIMENSIONS 


Intercooler Precooler 
0.40 0.40 
0.50 050 
0.60 0.60 
0.20 0.20 
0.05 0.05 


1.00 


TABLE 3-8 


HEAT EXCHANGER AREAS AND ASSOCIATED QUANTITIES 


Intercooler Precooler 


Regenerator 


Hr 0.2922 ‚3321 0.3321 


BL 0.07854 1047 0.1047 


0 


ur 


220 

. 00277 
.000491 
. 003649 


20379 


1774 


0.2688 


1 


30 


2220) 
‚00363 
. 000873 
. 006014 


. 0437 


27027 


45 


0.20 


0.00363 


0.000873 


0.006014 


0.0437 


0.2405 


6.717 





-- s e» ap 





TABLE 3-9 | 


LIST OF CONSTANTS FOR THE REGENERATOR, PRECOOLER AND INTERCOOLER 


Reference 
Quantity * Regenerator Precooler Intercooler Equation 
Ev. 2.407 0.1388 0.127 (3:25) 
Re 
RIV, 276.1 700 1165 (3-15c) 
Re 
۷ 446.5 3360 3360 (3-15c) 
h 
Y, 685.7 1110 1759 (3-4b) 
h | 
chy يآ‎ 1605.7 224,000 224,000 (3-4b) 
Er. 0.1575 0.1558 0.1566 (3-17) 
Er. 15 0.04098 0.04098 )3-17( 
TN 0.0667 0.0471 0.0471 (3-10a) 
KERE, VY O 9.73 x10 (3-13) 
Br 15.53 x 10 4 2010 x 10 ^ 2010 x 10 ^ ES, 
Le, 0.0805 0.02137 0.01694 (35125 
D^ /N 0.00490 0.00765 0.00765 (3-27) 
5 5 5 ok 
NV 4.288 x 10 2.02 x 10 1.2753 x 10° (3-28) 


lu 


* all units are as indicated in the list of 
symbols at the ‚beginning of the Chapter. 
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TABLE 3-10 

















TED HEAT EXCHANGERS 


Precooler Intercooler Units 





ft/sec 


Volume 
E 0.76 0.715 1.05 psi 


AP 12.03 4.87 4.17 1 


10.6 21.4 24.5 in: 





Description 


annulus area 
force coefficient 
chord 


specific heat at constant pressure 


Last stage isentropic enthalpy change 
ratio n ti it 
First stage 


diameter 
blade force 
reheat factor 


conversion factor, 32.17 


total enthalpy 
static enthalpy 
blade height 


conversion factor, 778 


. last stage actual enthalpy change 
ratio ES Y 7 E T 


mass flow 
rotational speed 
no. of stages 
pressure 


radius 


in. 


BTU/lb - °F 


ad 


LIO “2 
lb, sec 


BTU/1b 
BTU/1b 
ine 


ft = Ib 
BTU 





psia 


in. 








or 


= 


t" 6 dR 


Subscript 


Description 

entropy 

blade spacing 

absolute temperature 
peripheral velocity 
absolute velocity 
specific volume 

relative velocity 
absolute flow angle 
relative flow angle 
increment 

drag force/lift force 
efficiency 

density 

defined by equation (4-8) 
defined by equation (4-7) 
degree of reaction 


rotational speed 


Description 


compressor 
center, mean 
drag 

drag profile 
exit 


inlet, initial 
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Units 
BTU/lb - °F 
in. 

ºR 

ft/sec 
ft/sec 

ft /1b 
ft/sec 
degrees 


degrees 


rad/sec 











Subscript Description 


L Lift 

m axial 

n last stage 

R rotor, relative 

5 stator 

S stage, isentropic 

th theoretical 

u tangential 

x arbitrary stage 
General 


Before proceeding with a design study of the two compressors, the 
type of compressors must be selected. From Table 2-2 the energy impart- 
ed to the gas in the two compressors is 147.3 and 158.2 BTU/lb respect- 
ively. Analysis reveals that each compressor, if centrifugal, would 
consist of two or more stages in order to prevent excessive peripheral 
speeds. Rather than employ several centrifugal stages where efficien- 
cies are usually low due to losses occurring between stages, two axial 
Flow compressors have been selected. 

The flow channel associated with an axial flow compressor is 
illustrated in Figure 4-1 and the corresponding h-s diagram in Figure 


4-2. 
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Entropy 


Figure 4-2 Typical Enthalpy-Entropy Diagram for an Axial 
Flow Compressor 








From Figure 4-2, where it is assumed that the flow in the inducer 
and diffuser sections is steady and adiabatic, the following equations 


are obtained if geopotential energy differences are ignored: 





y? 
H 299 + Hh = constant (4-1) 
(Ah); 


(4-2) 
(A h s) Ul 


> 
IN 


ta 


(a hs) ne 


(4-3) 
(Ahlne 


IN 


a 2 
(^H), = (A و و‎ + xS T Ed au 


3 


2 


2 


G 
ro 

«o 
q 





| 
| Compressor Stage Performance 
| 


Before the compressor design can be carried out the type of stage 
must be selected. In this design only one type of stage was considered; 
namely, a non-symmetric stage where the absolute velocity entering the 
rotor is axial. Figures 4-3 and 4-4 show this type of stage and a 


corresponding h-s diagram. 
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Figure 4-3 Schematic Representation of an Axial Flow Compressor 
Stage and Corresponding Velocity Diagram 
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H, 
Entropy 
Figure 4-4 Typical h-s diagram for an axial flow compressor stage 


For the purpose of developing convenient working equations for 
calculating stage performance the velocity diagram of Figure 4-3 is 
changed to a dimensionless form by dividing all vectors by U to obtain 


Figure 4-5. It is assumed here that Ui = U, - 
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Figure 4-5 Dimensionless velocity diagram for axial flow 
compressor stage, 








In Figure 4-5 several new terms have been introduced. They are: 


C == Drag Force 
۴ Lift Force for rotor blade (4-5) 





Drag Force 
Lift Force for stator blade (4-5a) 


"| 





Wa, + Wa, 
ZU 


(4-6) 


vin 





EN 
Il 
4 
| 
ck 


== (4-7) 
T= Vaz = Va, E Wu, — Wu, (4-8) 
U U 
or 
_ = va for X, = 0 (4-8a) 


Expressions for stage efficiency are based on airfoil theory 
assuming incompressible flow and will not be derived here. Referring 
to Figure 4-5; if the product er E q is ignored with respect to En 


or e 


2p — U^$t | I-x-&d x -€ed| ار‎ 
is ^" 2 Eae ۱ T ri K 





س TE eses‏ هیا 


E = $ 1-2-6: , X-—Ée? | (4-10) 
Ts (R- p $t €s (1-29 dT Ce 
where — Cp (4-5) 
€ C. 








2 
e a 2020 + + 0.018 C. (Gan) 


S 
C, = == (4-12) 


NE 
C = L (2); (4-12) 


Equation (4-11) was obtained from a and is based on 
cascade tests and experimental single-stage and multiple-stage compres- 
sor data. These tests also show that the profile drag coefficient, 


C is primarily a function of the ratio s/c. This relation is plotted 


DP^ 
in Figure 4-6. (Howell) 

Figures 4-7 and 4-8 contain plots of additional test Fes co M 
showing a relation between fluid deflection, exit flow angle and ratio s/c. 
In addition to the stage efficiency, calculation of energy input 


per stage is necessary for determining the number of stages required. 


This can be expressed as: 


Ag, = 37 Ma 4) = 37 G 
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u ERU: 
For Vi = Va: 
E 
,و(و۵)‎ = port U 


The degree of reaction is defined as the ratio of actual enthalpy 


rise in the rotor to actual enthalpy rise in the stage. 


ha-h — hecho — (ama (4-15) 


h3-h, Ha- H, 5 (AH) 


Degree of reaction 





3, | 


If the relative flow in the rotor blade is steady and adiabatic: 


U 
H a E W "ED = constant (4-16) 
R h 29J ك‎ 
For U, = U, and Y = Wo: 

1 2 ma m, 

2 2 2 2 
h-h= Wi 7 Wa = Buc Wu (4-17) 
e | 2q J 2 9 J 
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From (4-8), (4-13), (4-15) and (4-17): 


Wa, + Wuz: 
2U 


Degree of reaction = 


x (4-6) 


Determining Number of Stages 


In a multistage compressor, the stage performance varies from 
stage to stage in some continuous manner. By evaluating the perform- 
ance of the first and last stages where the stream properties are known 


and assuming some reasonable relation between these and the other 





stages’it is possible to express U, T and Ms for any stage in terms 
of the two known stages. 


The relations assumed in the analysis are: 


Al 
0-۱ 
Ux — [oa (4-18) 
sf o 
a 
a A E e (4-19) 
x -1 


dx — e i (4-20) 
/ I 


From '4-l4a) it follows that: 


X-1 
UU. E: ABMS 
(ahy), oes ke (ah; (ahs), 7 





1 5 


۱ 


۲( ye 
(sho), Es ( * f)(AHj, = any pu yo" (4-22) 
X=2 


2 
where D = lahj, — (3 kn n (4-22a) 











o (A he), 5 ,لا‎ E N 
or 
n-2 
(1 +f) (AHs)in +) (A Hobin = n zx sper p B (4 -22b) 
(ahs), 075-۱ 
where n = number of stages required if for 
all stages (4h), = (Ah), 
Similarly: 
Ut (Ah) T 
ah = — T, = (ah) [46h a 
( he gy X ( ), H3 
4 
1 (Ah), a (A H), n ER 
X=] 
n-2 
(AH) : سا‎ | 
A hn == nu | + L a pe = (4-25) 
(Ah) ۱ 
۱ 2 
where = (Ahn — Un 0 (4-25a) 
(Ah), | t, 


3) 
and N = Number uf stages required if for 
all stazes (Ah), = (Ah), 
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Since there must be an integral number of stages the solution of 
equation (4-22b) for (Ah), generally will not be the same as the 
value obtained from the first stage caiculations. Thus a correction to 
the peripheral speed is usually necessary. Unless this correction is 
large, the effect on stage efficiency will be small. The total energy 
input can then be calculated with equation (4-25). 


The overall compressor efficiency is defined as: 


(A Hs) e (4-26) 


2 

where Cr Cp Ti a = (1-6a) 
2 2 2 

s (AH)ie = (AH), + ۷۰-۷۶ MM (4-42) 


N 29J cod 


The reheat factor, f, appears in equation (4-22b). Equations (1-8) 
and (1-9) give the reheat factors for an infinite number of stages for 
both expansion and compression processes. Referring to Figure 4-9 the 


reheat factor for hM stages of compression is obtained from: 


n 
Ah 
rf = و‎ (4 (4-27) 


he aie 
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Temperature, Enthalpy 





Entropy 


Figure 4-9 


For an infinite number of stages (4-27) becomes: 


Y dh: 
x = 
Ne ^ h. 


T (4-27a) 


Since dh = Tds + vdp for either a reversible or irreversible 


process and dh = Tds when P = constant: 


S, 
(Ah) = Tds = area under « constant pressure 
P. 


= Con si. 
mes 5 line on T-S diagram 
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| (Ah), = area between the stage inlet and 


= — 


exit constant pressure lines on the 
T-S diagram extending to a reference 


absolute temperature. 


3 3 
The difference between = su and = Celi is the small 


triangular areas shown in Figure 4-9. Assuming the constant pressure 
lines to be parallel and (Ah), constant for all stages, the area of 
i-f'-f consists of ۳ small triangles. 


Equations (4-27) and(4-27a) can be combined as follows: 


t Qo00 heohi > y GUD 
fo Mh- dhs 


x=] 


i 


area i-f-f' less n small triangles 
area i-f-f' 


or E 
+ 


oQ 
= n’-n = n-i 
n? n 
i n-li -2 
£ fee ^ (4-28) 


Compressor No. 1 


Having previously decided on the type of blading it is now neces- 
sary to determine an optimum combination of the many variables involved. 
First of all, a high stage efficiency is desired. To calculate this, 


practical values for (s/c) s (h/e), and (h/e), must be selected from 
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experience. The quantities, Cop? (s/c), and 4/8, can be obtained 
from test data of Figures 4-6, 4-7 and 4-8. Then for given values of 
e. the remainiug quantities in the dimensionless velocity diagram, 
Figure 4-5, and the stage efficiency can be calculated. 

The effect of a bad choice of (h/c), and (n/c)p is not great, 
since these quantities change the stage efficiency slightly. The 
‘values assumed for a preliminary investigation were 3.0 and 2.0 respect- 
ively. The choice of (s/c), is more critical. From a consideration of 
dimensions, especially blade heights, (s/c), of less than 1.0 proves to 
be impractical. Therefore, stage efficiencies were calculated for 
(s/c), of 1.0 and 1.5. Figure 4-10 shows the plot of fm vs. O, for 
these two ratios. It can be seen that the maximum efficiency in both 
cases will occur at some smaller angle than has been plotted. At the 
smaller angles the blade heights become smaller than practical design 
permits. In additíon, consideration must be given to Compressor No. 2 
and Turbine No. 1 since their rpm will be the same as Compressor No. l. 
For the same blading in Compressor No. 2 the smaller specific volume 
will result in even smaller blade heights. Too high a rpm will result 
in a large through flow velocity in Turbine No. 1 and consequently a 
low turbine efficiency. 

In conventional designs using air as the working fluid, the maxi- 
mum peripheral speed is usually determined from Mach effects because of 
the relatively low velocity of sound in air. The velocity of sound in 
helium is approximately three times as great so that allowable stress 
becomes the controlling factor. The materials that are presently being 
used for compressor rotors and blades will permit a U of about 


1200 £t/sec. 
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Figure 4-11 shows plots of Nvs. 4, for last stage blade heights 
of 1.0 and 1.4 inches. For calculation of the last stage blade heights 
it is assumed that the same blading is employed as in the first stage. 
From the plot it can be deduced that for given values of An and m 
the o sor having the smallest value of (s/c), has the lowest rpm. 
The compressors in this cycle require high rpm if small blade heights 
and large diameters are to be avoided. High Ere rpm results in a 
blade height to mean diameter ratio which is too large, so some compro- 
mise has to be made, Preliminary calculations indicate that 7000 rpm 
is satisfactory for Turbine No. 1 and the two compressors. Another 
factor to be considered is the number of stages required for different 


values of (s/c),. 


TABLE 4-1 


EFFECT OF (s/c)p ON COMPRESSOR BLADING 


(s/c)p 1.0 1-5 Units 
N 7000 7000 rpm 

E > 02.2 64.4 degrees 
h. 1.4 1,4 ” inches 
Ns 2952 .873 2 

n 7 or 8 10 stages 
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From the comparison in Table 4-1 the blading having a space to 
chord ratio of 1.0 appears to be the most favorable and is used in 
subsequent analysis. 

The results of preliminary calculations indicate that an (h/c), 
or 3 results in a very small chord so it was changed to 2.0. 


In order to make all rotor disks of the same dimensions, the blade 





root diameter will be constant. The first and last stage data that 


result are Listed in Table 4-2. 


TABLE 4-2 


DATA FOR FIRST AND LAST STAGES OF COMPRESSOR NO. I 


First Stage Last Stage Units 

X, 0 0 degrees 
(Ya 1.0 1.0 € 
GAR 2.0 1.63% a: 
(hk). 20 1.63* coos 

6, 62.0 62.0 degrees 
0 70.6 70.6 " 

d 0.3522 0.3522 LT 

E 0.3376 0.3376 ose 

2 0.8312 0.8312 ies. 
x, 43.8 43.8 degrees 
(39). 0.48 0.48 e 

N: 0.8812 0.8755 ue 

Er 0.0482 0.0513 nn. 

Es 0.0438 0.0440 mæ 

* h/c for last stage = h/c for first stage x i 
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Since there must be an integral number of stages the calculation 
of work imput per stage, dimensions, and velocities requires a series 
of approximations. Table 4-3 containg, first approximations of U, و۷‎ 
A , h and d, for a fixed compressor speed of 7000 rpm. Before 
proceeding with the analysis, some additional information is necessary; 
namely, the inlet and exit pipe velocities and the inducer and diffuser 
efficiencies. The velocities were selected after taking into consider- 
ation the size of pipe required and the Uae are typical values 
E conventional designs. These are indicated in Table 4-4 along with 
calculated results. 


The following equations, supplementary to those listed previously, 


are also applicable: 





| 


A IX 144 — mhd in? (4-29) 





rpm (4-30) 


TABLE 4-3 


FIRST APPROXIMATION DATA FOR FIRST AND LAST STAGES OF COMPRESSOR NO. 1 


‘First Stage Last Stage Units 
N 7000 7000 rpm 
U* 1200 1190 ft/sec 
VQ? 423 420 ft/sec 
A * 213 172 in 
h* 1.72 1.40 in 
E 99 92 39.00 in 


سا 


* First approximation 
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nº 


(ahs) 

(AR), 

(AH), 
Tun 
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* First approximation 


BTU/1b 


t 








Therefore, some changes. 


Corrected results for the first and last stages are indicated in 


Table 4-5. 








TABLE 4-5 


FINAL DATA FOR FIRST AND LAST STAGES OF COMPRESSOR NO. 1 


First Stage Last Stage Units 
U 1235 1225 ft/sec 
N 7000 7000 rpm 
Vi 436 432 ft/sec 
Vo 604 399 ? 
8 1310 1300 E 
W, 928 920 E 
A 207 167 in? 
h 1.63 1.33 in 
d, 38.87 38.87 C 
d. 40.50 40.20 i 
d. 42.13 41.53 P 
c 0.82 0.82 9 
Sp 0.82 0.82 " 
So 0.39 0.39 S 


Compressor No. 2 


The same blading is used in Compressor No. 2 as was used in Com- 
pressor No. 1 and the rpm is fixed at 7000. Figure 4-12 shows the 
variation of A, and n as a function of U. A mean peripheral speed of 
1050 ft/sec was selected in order to obtain a satisfactory last stage 
blade height. The root diameter is constant as in Compressor No. 1. 

Table 4-2 is applicable with the following exceptions for the 


last stage: 
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Calculated results are listed in Tables 4-6, 4-7 and 4-8. 


TABLE 4-6 . 


FIRST APPROXIMATION DATA FOR FIRST AND LAST STAGES OF COMPRESSOR NO. 2 





First Stage Last Stage Units 
1050 1037 ft/sec 
7000 7000 rpm 
9 370 365 ft/sec 
A 164 134 in? 
h* 1252 1.26 in. 
0.5 34.38 33.90 


elo 


' First Approximations 
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TABLE 4-7 


ADDITIONAL DATA FOR COMPRESSOR NO, 


128.83 
13.12 
2599 
1.0157 
10 
0.9793 
0.9849 
9.9005 
9.9189 
13.20 
14.98 
148.58 
.867 


. 840 


ft/sec 


BTU/1b 


n 


stages 


BTU/1b 
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TABLE 4-8 


FINAL DATA FOR FIRST AND LAST STAGES OF COMPRESSOR NO. 2 


First Stage Last Stage Units 
1050 1037 ft/sec 
7000 7000 rpm 

370 365 ft/sec 
513 507 H 
1115 1100 E 
788 778 H 
164 134 in^ 
12752 1.26 in. 
32.86 32.64 N 
34,38 33.90 Ç 
36.90 35.16 H 
0.76 0.76 £ 
0.76 0.76 ۱ 
0.37 0.37 
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CHAPTER V 








TURBINES 
Symbols used in this chapter: 
Symbo 1 Description Units 
d diameter ing 
f reheat factor Bern 
g conversion factor, 32.17 a e 2 
lb, sec 
H total enthalpy BTU/1b 
h blade height in. 
or 

enthalpy BTU/1b 
J conversion factor, 778 a 
K coefficient defined in equation (5-10b) --- 
m mass flow Lbs/sec 
N rotational speed rpm 
n number of stages --- 
P absolute pressure psia 
r radius in. 
T absolute temperature °R 
U peripheral speed ft/sec 
V absolute flow velocity ft/sec 
V Specific volume ft MIB 
W relative flow velocity ft/sec 
OL absolute flow angle degrees 
8 relative flow angle degrees 
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LD 
KT 


ol Description 


increment 
efficiency 

defined by equation 
defined by equation 
coefficient defined 
defined by equation 
degree of reaction 


coefficient defined 
and (5-13) 


c  —— 


rotational speed 


Subscript Description 

B blade 

E center, mean 

e exit 

i inlet, initial 

m axial 

n last stage 

0 outer 

R relative 

S isentropic 

T turbine 

u tangentíal 
General 


(5-16) 
(5-15) 
by equation (5-7) 


(5-20) 


by equations (5-12) 


Units 


rad/sec 


This chapter consists of a design study of the two turbines in the 


cycle. As in Chapter IV for the compressors, a one-dimensional flow 
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analysis is performed to determine turbine diameters, blade heights, 
number of stages and blade shapes at a mean streamline. The flow is 
assumed to be steady and adiabatic in both moving and stationary 
channels so that the total enthalpy, H, in the stationary channels and 
the total relative enthalpy, Hp» in the moving channels are constant 
along the streamlines. Analysis reveals that if a centripetal turbine 
is employed an excessive peripheral speed will be required to produce 
the necessary work output. Therefore, the entire analysis that follows 


is based upon an axial flow turbine. 


Stator 


Inlet Pipe Exit Pipe 
lo 
L 
f 
Rotor 
Figure 5-1 Schematic Representation of an Axial Flow Turbine 
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where hi = 





where H = 


| 
m 

| 
AJA 
5 


AHg E 


(&Hs) — (A is Le A 


(A o | 
and nr — Il 





ate and accelerate the flow respectively, since, 


Figure 5-2, the following relations can be written: 





Figure 5-1 is a sketch of an axial flow turbine showing stream 
locations i, 1, n and e. The inlet and exit sections serve to deceler- 
as we shall see later, 
the pipe velocities for this cycle are necessarily larger than the 


through flow velocity. Referring to the enthalpy - entropy diagram in 


(5-la) 


(5-1b) 


(5-2) 


(5-2a) 


(5-3) 


(5-4) 





Enthalpy 


3E 
de 

JE 
© 





Figure 5-2 





Entropy 


Typical Enthalpy-Entropy Diagram for an Axial Flow Turbine 
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Stage Efficiency 


Figures 5-4 and 5-5 are enthalpy - entropy diagrams illustrating 
absolute flow quantities for the stage and relative flow quantities for 
the rotor blade respectively. 


From Figure 5-4: 


۲ 
Ah) = hp-hz? (5-5) 
(ahs) = Ms ha? (5-6) 
(ah) = h,x -ha» | (5-7) 


Assuming zero reheat (f = 0): 


(Ah.)" "M (A he) (EN 


d EE (Ah)vz (5-8) 


(hs)? 
9) 
x? = (Ahs) = isentropic degree (5-9) 
ARV of reaction 
— Ah 
na 28, = u ,x’ (5-10) 


(ahd FOR 
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(a) 


Figure 5-3 


HH, 


Enthalpy 


Figure 5-4 


2 
E. 





Stator 


T€ 





fc 
f 
(b) 


Velocity Diagram for an Axial Flow Turbine Stage 








Va 
R r 29 ah, 
+ (ah) 


AD 
25 


Entropy 


Typical Enthalpy-Entropy Diagram for an Axial 
Flow Turbine Stage 


88 





In equation (5-10) the kinetic energy of the fluid at the exit of 


the stage is not considered as useful energy. Actually, most of this 


energy is available for the succeeding stage and another stage efficien- 


cy is defined: 


nã = aay (5-10a) 
2 
where (A he). ze (sh) + K Ma (5-10b) 


29J 


and K = percentage of the kinetic energy 
leaving a stage that is available 
for useful work in the succeeding 


stage 
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Enthalpy 





Entropy 


Typical Enthalpy - Entropy Diagram for Relative Flow 


Figure 5-5 
in a Turbine Rotor Blade 


From Figures 5-4 and 5-5: 





= = tant 5-11 
Ar, Hr, constan ( ) 
E w^ | U* ۱ 
where Hke = h+ 293 293 (5-lla) 
Y’=|ı- Reeth (5-12) 
2 WE 
3%93 
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Y = Ws (5-13) 


W397) 
2 2 2 2 
Ve 2 Do ya We +U3-U, 
Na س‎ A A A 2 2 


For zero reheat (f = 0): 


) 
2 7 2 2 £ 
۱ 22 و‎ h 
Us 3 ê ie 
2 EE 2 
It is advantageous for subsequent analysis to introduce the 


following definitions: 


2 
2 GE 


Y Lc ۳ 
(ahs) (293) “a 





L= E (5-16) 
2 


Vm 
Vm 


3 (5-17) 





H 


2 


From the trigonometry of Figure 5-3, the definitions above and 


equation (5-l4a), several new equations can be derived: 
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qj - Gil «even - 


2 ۵ (۶ ۲۷۱-۵ sink, + سر‎ 


۷ 

Lr vs (5-18) 
A Y A 

ie — e v (6 VI-X? Sind, + AY X- x») (5-19) 


Y= X" CLLKT —- e cos* oc (1- X) - 


2 
20 4 Y NIT 5۱ ولا‎ + y^ (* E ) (5-20) 


A c (5-21) 
tan 8, ٩5 7 cos A, 


The work output per stage can also be expressed in the familiar 


form: 


۸۳ == o LA RV.) (5-22) 
TE R = e 

pU | 
Ah, — 97 M. Vas) (5-22a) 





If x, = U 


V 
Aha = 33 | (5-22b) 


E 





When all stages are similar and produce equal work: 
Nas DPA, (5-22c) 


From the definition of reheat factor and equation (5-10b): 


2 


) 
(AH), mm n(ah) ^ jh ahe 7. 


i qe 20 ل‎ +۴ 


where f= fÉ A (4-28) 


Two additional relations necessary for turbine design are: 


— 720 U - 
d. EA IT N un. (5-24) 


Den TR in? (5-25) 


m 


Selection of Blading for Turbine No. 1 


The type of blading to be selected must permit a good efficiency 
as well as reasonable dimensions and number of stages for the condi- 
tions imposed. The use of low efficiency Curtis impulse stages would 
not be satisfactory despite the fewer stages required. Generally, it 


can be expected that reaction stages will have the highest efficiency 
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and require the greatest number of stages. In order to obtain a high 
efficiency and prevent too many stages a compromise between impulse and 
reaction stages may be warranted. Two general types of stages were 
investigated; namely, (1) symmetric or 50 per cent reaction and (2) 
modified Rateau. The modified Rateau stage differs from the pure 
impulse stage in that a small amount of reaction is introduced by 
specifying an axial exit velocity. 

Comparisons of the two types have been based on conditions exist- 
ing at Turbine No. 1, the rotational speed of which is 7000 rpm. The 
maximum peripheral speed, U, is limited by stress conditions. At 
temperatures of 1350°F a peripheral speed of about 800 ft/sec seems 
possible with materials currently available. Losses in the stator and 
rotor, expressed by the definitions of d ۲ Yo: and Y'a are esti- 
mates based on conventional gas turbine performance tests. For a 
single stage the differences between U, and U, can be considered negli- 
gible. In addition, the tifrough-flow velocity, Vo? is assumed to be 
constant. 

For each value of ue there is an efficiency curve which can be 


plotted as a function of — . The choice of ve depends to a large 
2 


extent on the ratio of blade height to mean diameter. From Table 2-2: 


E SS ental: 
m, — 182.4 Ib/sec 


VE 2019 ft/lb 
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From equations (5-24) and (5-25): 


d, = 2765 in. 
Vinh, = 15176 f#t-ın 
Sec 


K = 300 ft/sec results in a first stage — ratio of 0,13. Velo- 
c 


cities less than 300 ft/sec result in higher efficiencies but the blade 


heights become too large, especially at the last stage. 





(a) Modified Rateau Stage 





(b) 50% Reaction Stage 


Figure 5-6 
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Using the quantities assumed below, the number of stages, n, and 
stage efficiency, ne were calculated for the two types of stages 


indicated in Figure 5-6. 


N = 7000 rpm p = 0.95 
U = 800 ft/sec Y = 0.95 
Mz 1.0 Y = 0.92 
A = 1.0 K = 0.90 


Plots of n. and n vs. - are shown in Figures 5-7 and 5-8. 
2 
It is apparent that there is little difference in maximum efficiency 
and that two less stages are required for the modified Rateau stages. 


For this reason the modified Rateau stage was selected. At the condi- 


tion of maximum efficiency: 


X- = 0.7 
2 
Ahg= 35.3 BTU/lb 
* 
s = 0.854 
5 
x. = 099) 


n = 8,7 theoretical stages 
An adjustment must be made in the work per stage to increase the 
number of stages from 8.7 to 9. For a nine stage turbine having equal 


work per stage: 


(ah) — SHe — 340 BTU/Ib 
BJ actual 


9 
UR O ju. l 
Since ۰ = “oF = G3 (e) (5-22b) 
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i 
1 
L 


1 
T 
s 











UN = Va (Mas /U) x (Aha) actual 


(Va, / D ۵ 


Ya, 
To use the same blading, ۳۳ ۳ constant. 


Io 


( ) 


(DA PORT — EP 800 = 785 ft/sec 


, 


Recalculation of velocities and dimensions assuming that the 


blade root diameter, d is constant along the turbine, produces the 


1" 
results contained in Table 5-1. 
In calculating the turbine efficiency, 7 po inlet and exit effects 


and reheat factor must be known. The reheat factor can be calculated 


by the relation: 





f= fo 2 (4-28) 


where foo is obtained from equation (1-9). 
To obtain (4 Hs); e , the values of Vo Vos /FL and H. must 


be known. The pipe velocities that are assumed here will be covered 


iu more detail in Chapter VII. 


Vi = 600 ft/sec 


e = 500 ft/sec 
Ai z 0.85 

i 
n = 0.75 

e 
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BTU/1b 


BTU/1b 


PEU. | 
Recalculation of the exit pressure and specific volume is not 
¥ 
2ssary since the calculated efficiency is only slightly higher than 


> assumed value. A summary of the first and last stage data for 


bine No. 1 is listed in Table 5-1. 








۳ E are LO 19 


FIRST AND LAST STAGE DATA FOR TURBINE NO, 1 





First Stage Last Stage Units 
1121 1162 ft/sec 
294 305 " 
74. 74.78 degrees 
0 " 
ft/sec 
W, 847 878 " 5 
6, 45.29 45.29 degrees 
ei -67.98 -67.98 " 
U 785 613 ft/sec 
di 21.89 21.89 Ln 
d. 21215 28.10 M l 
d. 32.41 34.31 " 
h 5.26 | 6:21 " 
h/d, .194 221 --- ] 


Additional Turbine Data: 


= 
u 


7000 rpm 
x 

Ds = 0.854 

Ir = 0.858 


9 stages 


- 
H 


101 








Turbine No. 2 

Unlike Turbine No. 1, this turbine is not restricted to a specific 
design rpm, thus the speed of this turbine is an additional variable. 
Since the conditions for this turbine are not very much different from 
those of Turbine No. 1, a similar type blading is used for simplicity. 
Thus, U, Va? M AE 2 : Y >: y^» and K remain the same. From 
Table 2-2: 

AHg = 193.7 BTU/lb 


182.4 lbs/sec 


De 
ü 


VQ = 6.692 £c? /1b 


720 U 
S — T L 5-24 





h = MANxI44 _ 0.001017 N inches 
Vm, U, x 720 


For N = 5500 rpm: 


hy = 5.59 inches 


a er 


"1/8, ROMS | 


1 


The theoretical number of stages required is 5.5. Using the same 
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method as employed for Turbine No. 1: 
h = 6 stages 
V. = 500 ft/sec 


V = 625 ft/sec 


e 
i = 0.85 
lee = 0.75 


LL = 1.00825 
(A Hs) „= 226.5 — BTU/1b 
4 


(A Hs), = 231.8 ۰ ۵ 3 
na = 0.835 


Since this efficiency is 3.5% larger than assumed in the prelimin- 
ary analysis, a recalculation of exit pressure and specific volume is 


necessary prior to calculating blade heights and diameters at the last 


stage. 
W 
$) co = (AHI Le 
Zell] (1-7a) 
R Ty Cp 


e s — 4576 psia 
V = 8.313 ۸ 


A summary of first and last stage data for Turbine No. 2 is listed 


in Table 5-2. 
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First Stage Units 





1095 | ft/sec 
287 i É 


74.18 degrees 


degrees 


ts 


ft/sec 





Additional Turbine Data: 


N = 5500 rpm 

NÝ = "0.854 
N = 0.834 

n = 6 stages 
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CHAPTER VI 


REACTOR 
Symbols used in this chapter: 
symbol Description 
A total heat transfer area 
a total coolant flow area 
C fission rate per unit of power 
E specific heat at constant pressure 
d diameter of coolant channels 
f thermal utilization factor 
g conversion factor, 32.17 
H height of reactor core 
h convection heat transfer coefficient 
3 size-shape factor 
k' thermal conductivity 
kæ multiplication constant for reactor 
of infinite size 
Koff multiplication constant for reactor 
of finite size 
L thermal diffusion length 
L thermal diffusion length for pure material 
M mass 
m mass rate of flow 
N number of eed eae 


or 
number of coolant channels 
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Units 


ft? 


er 


fissions/watt - sec 
BIU/1b °F 


in. 


lb = ft 
m 


Ib, = eer 


cm or ft 


2 


BTU/hr - ft" - °F 


-2 
cm 


2 


BTU/hr - £t” - *F/ft 


cm 
cm 
grams or Ib 


lb/sec 


< 


< 


BP “4 مهم يت‎ n kx A = 





Description 


power 

pressure 

resonance escape probability 
rate of heat transfer 

radius of reactor core 
temperature 

reactor volume 

specific volume 

coolant velocity 

microscopic cross-section 
absolute viscosity 

fast fission factor 

slowing down area for fast cons 
neutrons born per absorption 
flux 

friction factor 


increment of 


subscript 

a absorption 

av average 

E graphite 

f fission 

He helium 

m mixture of graphite and uranium 
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Units 
Watts 


psia 


BTU/hr 
Cl or CE 
°F 

one or ft? 
£t? /1b 
ft/sec 


cm 


lb/hr - ft 


2 
neutrons/cm - sec 








Subscript Description 


u uranium 


w wall 


General 
Although the writer has had límited training and no practical 
experience in the field of reactor technology, an attempt was made in 
this chapter to determine the approximate dimensions of the core of a 
helium-cooled enriched reactor, for the primary purpose of establish- 
ing the pressure drop of the e 
Critical Mass 
In determining the core dimensions and the amount of fuel required 
to produce a critical mass the following assumptions have been made: 
(1) The core consists of a homogeneous mixture of graphite and 
235" 
(2) The core shape is a circular cylinder with a height, H, 
equal to its diameter, 2R. 
(3) The coolant channels occupy 1/3 of the total volume and are 
uniformly distributed throughout the core. 
(4) The power rating of the reactor is based on a cycle thermal 
efficiency of 0,329, 
(5) An average flux at full power is of the order 0 neutrone a 
- sec. 
(6) The effect of reflector and cooling channel liners is not 
included. 
(7) No allowance is made for fuel depletion, poison and tempera- 


ture effects on reactivity. 
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eee 
- C 
Ke € K 


72 ag 


Kem =| = 


BN — rr 


= چ‎ + GY 





A 


(É = ES (1— f) TL 


| 
f — NG (6-5) 
Nu (alu 


To calculate these quantities the following nuclear data is known: 
(6 a), =z 0.0045 barns (1 barn a jo ^^ ema) 
(Ca), = 650 barns 


p = 2500 em 
w 
p ^ 1.0 
7 
Rene 
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\ 


Preliminary calculations indicate that for the power and flux 





requirements: 
we 
mr / 

N حر‎ 14,400 

u 
Therefore: 

f = 0.909 

kæ 1.91 

E 


By trial and error calculations using equation (6-1): 
v = 0,0013 


H = 159 cm = 5.23 ft 


To allow for the effect of the cooling channels and considering 





that they are distributed uniformly over the entire core, the following 
approximations have been made: 
(1) L? increases to 9/4 of its original value 
(2) T increases to 9/4 of its original value 
Referring to equation (6-1) the net effect of Ke is to decrease it 
to 4/9 of its original value or: 


K? = 0.000578 
E 02) cm$ Ew 65 ft 
The ratio of graphite to uranium by mass is: 
DOZEN EN 3 
Mu 235 Na 
Since it can be expected that the density of the graphite will be 


Z ۱ 
somewhat less than its pure state (1.65 g/cm ), it is assumed to be 
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equal to 1.55 Em. 





Then: 
x H? 6 3 3 
we æ 2/3 7 TO 10 CM E 7253 ft 
H. = 7 10º xo ce 11.15 x 10º grams = 24,600 lbs 
6 
۱۷ ۱ 
M, = 736 = 15150 grams = 33,4 lbs 
1.55 x 6.023 x 10% 22 3 
c 12 
22 
7.77 x 10 18 , 3 
N, = 14,400 = 5,4 x 10 nuclei/cm 


The power rating of the reactor, based on a 50,000 HP output at a 


thermal efficiency of 0.329, is: 


50,000 


pus سح‎ SOOO 
P 0.329 ó 
ES SS SS Walls 
The average flux is: 
Zu c = neutrons 
Po EV e 
i Na 7 Vo, cm - sec 
where Go EEE RX 10° fissions/watt - sec 
O E barns 
6 3 
MEE 6/75 x 10 em 
or 
d = 1 6 x ee neutrons 
cm ~ sec 
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The burn-up rate is: 


No. of vae‏ تخت بر ی ی 


U 
day fission absorption 


=235 


Ja 235 
c x 3600) P A So o 
( (G), 6.023 x 10º? 


136 grams/day at full power 


Cooling System 


Pd 


Although the dimensions of the core and the amount of fuel and 
moderator mixture will increase when allowance is made for the effects 
of fuel depletion, poison and temperature, it is assumed here that this 
will have little effect on the characteristics of the cooling system. 
From the helium properties listed in Table 2-2 for the inlet and exit 
of the reactor and the dimensions of the core: 


Mae eV = és qt 


e 


7 Me — 6.15 ft 
H 
TA = na = 49 ft/sec 


For an average wall temperature of 1500°F the heat transferred to 


the coolant is: 


O A(t, - es = Cpm (t.e tin) x 3600 (6-8) 
BTU/hr 
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where ۳ E = 387.8 °F 
wW av 





t Sb = 475.5 °F 
out ın 


m = 182.4 lb/sec 


1.24 BTU/Ib - °F 


c 
0 
H 


or h 





L 0 - 'F‏ ا 
A‏ 


_ 4 . 2960 
N = ۲ ۰ re 


Td — 6090 frt 
DONE. ONU Hom == tê 
le dl 


LXI NEmUhr .-ft'-*F (9 


Another expression for h as a function of d is obtained from 


the turbulent flow heat transfer equation: 


0.8 | 
Nees O 02 3 Ke me (3-16) 


in which the following gas properties are evaluated for an average 
coolant temperature of 1119°F; 


"Ecc esce is 


C = 1.24 BTU/lb - °F 


2 


M = 0.10 lb/hr = EC 


0.18 BTU/hr - ft - °F/ft 


n 
H 
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- °F 


b = un BTU/hr - ft 


2 


From equations (6-9) and (6-10): 


1.08 inches 


2 


178 BTU/hr = ft - °F 


2540 channels 


2 Pipe Friction 





o 
CEE = 00225 (3-14) 
where Ro = 36,500 


1 O os: 
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TABLE 6-1 


SUMMARY OF THE REACTOR CORE DATA 


H 7.85 ft 
M, 15.15 Kg 
H. 11,150 Kg 
۱ 
۷ 379 £c? 
W 49.0 ft/sec 
d 1.08 inches 
N 2,540 channels 
p 0.117 psi (Pipe friction only) 
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Symbol 
d 


< 


A 
n 
Be 
$ 
۳ 


Subscripts 


i 


O 


th 


CHAPTER VII 


INSTALLATION 


Symbols used in this chapter: 


Description 


diameter 


conversion factor, 32.17 


mass rate of flow 
pressure 

thickness | 
velocity 

specific volume 
increment 
efficiency 

stress 

defined on page 118 


defined on page 118 


inside, initial 
outside 
thermal 


wall 
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Units 
in 


ho = ft 
2 


ior ^ sec 
1b/sec 
psia 

in 

ft/sec 


LIO 








General 

In the previous chapters each of the major components of the cycle 

has been analyzed. The following data were obtained: 

(1) Basic internal dimensions of the turbo-machines such as 
diameters and blade heights. 

(2) Turbo-machine characteristics N as stage and overall 
efficiencies, number of stages, and velocity diagrams at the 
mean radius for the first and last stages. 

(3) Heat exchanger dimensions including dimensions and arrange- 
ment of tubes, shell inside diameter, and total internal 
volume. 


(4) Heat exchanger pressure losses due to pipe friction. 


In order to combine these components in an installation, addition- 
al information is necessary. Such items as piping, casings, inlet and 
exit connections, bearings, couplings, auxiliary equipment, control 
apparatus, and expansion joints, must be designed. It is beyond the 
scope of this investigation to cover all of these items, so only those 
which affect the size of the plant and pressure drops will be consider- 
ed in this chapter. 

Before considering the individual components there follows a 


review of some of the relations used in this chapter. 


‘> 
Effect of Changes of Specific Volume between Inlet and Exit of 


Heat Exchangers 


In the calculations for pipe friction using equation (3-13) a mean 


specific volume was used. Given the channel below: 








ER‏ سے ر 


Figure 7-1 Channel of Constant Flow Area 


A 
the pressure change between (1) and (2) due to the change in specific 


volume can be approximated by the equation: 


e 
E ~ Mn 5۱ 7-1 
gm SL qua, Po SSR 


Pressure Drop for Flow Acceleration 


ARNS NA 


Figure 7-2 Channels of Decreasing Flow Area 


2 We 
2 | ¢ \2av x 144 
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where Q is the efficiency of conversion from static 


“pressure to velocity head. 


Pressure Increase for Flow Deceleration Due to a Sudden Enlargement 


\ 





\ 
EF E 
e | 
a 
\ 
Figure 7-3 Channel with an Abrupt Increase in Flow Area 
From the momentum theorem, 
2 
۳ ۲ 2 ۱ 
R-P = 4241 psi (7-3) 


| gv x 144 


Pressure Increase in a Diffuser 


Y BRR Soa 


TR 


Figure 7-4 Diffuser 


2 2 
DEM yvy) pst (7-4) 
i O + 


4 


where Y is the efficiency of conversion from velocity 
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head to static pressure, 


Pipe Dimensions 


The internal pipe diameters and flow velocities are related as 


follows: 
L 
2 
— = (EY | x 12 inches (7-5) 


The wall thickness, t , depends on the allowable circumferential 


stress, f 


(^ P) d; 


inches (7-6) 
L0 








Piping 
In analyzing the turbines and compressors inlet and exit pipe 
velocities were specified in order that overall efficiencies could be 
determined. The internal pipe diameters corresponding to these velo- 
cities can be calculated.from equation (7-5). It is desirable to use 
standard pipe sizes although these may differ slightly from the calcu- 
lated diameters. Using the pipe property cnart and allowable stress 
data from Re o Table 7-1 was compiled with the use of equa- 
tions (7-5) and (7-6). 
The actual outside diameters of some of the installed piping will 
be greater than those listed by the amount of insulation required. For 
| locations four through eight (Figure 2-1) this will be in the order of 


six to eight inches. For the other locations none will be required 


unless some is desired to limit the space temperature. 
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Regenerator and Precooler 

It is convenient to discuss these two units together since they are 
in series for cooling the hot gas between the power turbine discharge 
and the inlet to Compressor No. l. Figure 7-5 illustrates the two units 
showing approximate dimensions. For an allowable stress of 10,000 psi 
for the regenerator shell, its thickness is two inches. An allowable 
stress of 12,000 psi for the precooler shell requires a thickness of 
one inch. For purposes of calculating pressure changes, schematics of 
the flow passages are illustrated in Figures 7-6 and 7-7. 

Referring to Figure 7-6, the flow has been broken down into the 
following parts: 

(a) Deceleration by means of a diffuser of 2 foot length between 

(1) and (2), with corresponding pressure recovery. 

(b) Deceleration to zero velocity at (3) with no pressure recovery. 

(c) Acceleration to E 

(d) Pipe friction effects between (4) and (5). 

(e) Deceleration from Vs to zero velocity at (6) with no pressure 

recovery. 
(£f) Acceleration to V 
(g) - (k) For remainder of the channel between (7) and (12) the 


effects are similar to those outlined in items (b) - (f) above. 


(1) Effects of changes in specific volume between inlet and exit. 


A summary of the velocities and pressure changes for the locations 
specified in Figure 7-6 is listed in Table 7-2. Where applicable, 


W - 0.85 and Ọ = 0.80. 
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’ 


Figure 7-6 Schematic Representation of Flow through the 
Regenerator and Precooler 





Tube Sheets 





Figure 7-7 Schematic Representation of Flow through the 
Regenerator (H.P. Side) 
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TABLE 7-2 
ica mde 15 


ARY OP PRESSURE LOSSES IN THE REGENERATOR AND PRECOOLER 













on (x) _ 
2 7-6) ver RE) V( ft/sec) E D) Equation 


۳ X nn... 





8.413 


458 


0 


il 


12 E 400 -6.30 (7-2) 
Compressibility Effect 0.12 (7-1) 
A "total ^ po 
EH = 0.0158 
= 
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TABLE 7-3 


SUMMARY OF PRESSURE LOSSES IN THE REGENERATOR (H. P. SIDE) 


Location (x) 


(Figure 7-7) CE /16) V(ft/sec) E ۱ e) Equation 
1 1.883 450 === ۳ 
2 y 323 1822 (7-4) 
3 x 0 0 (7-3) 
4 66 -0.31 (7-2) 
5 35677 132 -12.03 (3-13) 
6 2 0 0 (7-3) 
T " 600 -13.20 2) 
Compressibility Effect -1.09 (7-1) 
A P = -19 43 
Ar 
E a 0.0194 


Kas 


Considering the conditions inside the tubes of the regenerator, 
the schematic of Figure 7-7 indicates that items (a) - (f) on page 121 
apply. The results are listed in Table 7-3. Where applicable, 


WY = 0.85 and @ = 0.80. 
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Figure 7-8 Intercooler 





Figure 7-9 Schematic Representation of Flow through the Intercooler 
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Intercooler 

The intercooler is similar to the precooler and the flow is broken 
down in the. same manner. The dimensions indicated in Figure 7-8 are 
approximate and based on calculations in Chapter III. For an allowable 


stress of 14,000 psi, the shell thickness is one inch. 


TABLE 7-4 


SUMMARY OF PRESSURE LOSSES IN THE INTERCOOLER 


Location (x) 








(Figure 7-9) v(£t?/1b) V(ft/sec) Ex ۱ ar) Equation 
1 2.768 400 Sey a 
2 G 294 22465 (7-4) 
3 Ñ 0 0 (7-3) 
4 = 82 -0.32 (7-2) 
> 2.309 68 -1.05 (3-13) 
6 0 0 (7-3) 
7 B 300 -5.26 (7-2) 
Compressibility Effect 0.09 (7-1) 
AP = -4.11 
> Û 0062 


" 
ds 
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Reactor 

As illustrated in Figures 7-10 and 7-11, the flow through the 
reactor is similar to that through the heat exchangers. The overall 
dimensions are only estimates, since no calculations were made concern- 


ing the thickness of shielding and reflector required. 


TABLE 7-5 


SUMMARY OF PRESSURE DROPS IN THE REACTOR 


Location (x) 2:15) V(ft/sec) i ۱ P,.¡(Psi) Equation 
1 3.677 600 weer uu 
2 > 434 4,29 (7-4) 
3 D [578 0.34 (7-3) 
4 1 41.5 -0.06 (1-2) 
5 5.019 Sov» 0012 (3-13) 
6 Vi 18.9 05083 (7-3) 
7 É 600 -9.69 (7-2) 
Compressibility Effect -0.04 (7-1) 
RR. -5,25 
a = 0.0054 
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Schematic 


7-11 


Flow through the Reactor 





Tables 7-6 and 7-7 compare the calculated machine efficiencies 
and pressure losses with those assumed in Chapter II. Another cycle 
analysis using these calculated values results in the following new 


performance indices: 


162.5 lb/sec (For 50,000 HP output) 


m 


7] th = 0.361 


helium rate 11.7 1b/HP - hr 


i 


work ratio = 0.424 


TABLE 7-6 


COMPARISON OF ASSUMED AND CALCULATED VALUES OF AY /P; 
Unit ۵ (assumed ) A P/P; (calculated) 


Regenerator 


High pressure 


side 0.015 0.0194 
Low pressure 
side 0.010 
0.0158 
Precooler 0.0075 
Reactor 0.010 0.0054 
Intercooler 0.005 0.0062 
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TABLE 7-7 


COMPARISON OF ASSUMED AND CALCULATED VALUES OF MACHINE EFFICIENCIES 


n (assumed) n (calculated) 





Compressor No. 1 0.84 0.867 
Compressor No. 2 0.840 
Turbine No. 1 0.85 0.858 
Turbine No. 2 0.80 0.835 
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CHAPTER VIII 


STRESSES 


Symbols used in this chapter: 


Description 


area 

width of turbine disk 

chord 

specific heat at constant pressure 
force 


conversion factor, 32.17 


blade height 

enthalpy 

moment of inertia 
conversion factor, 778 
coefficient defined in equation (5-10b) 
bending moment 

mass 

mass rate of flow 
rotational speed 
pressure 

centrifugal Loading 
radius 

Spacing 


absolute temperature 
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Units 


in 
in 
in 
BTU/lb - °F 


lb. 
DE Et 


lb 2 


f sec 
in 
BTU/1b 
in’ 


ft - 1b/BTU 


in - lb 
1b 
lb/sec 
rpm 
psia 
Lb/in 
in 

in 


°R 








N 


R 


` 


o XT 
O 4 


ax do bE 1.0 e 


subscript 


eo 


CG 


Description Units 
temperature °F 

peripheral speed ft/sec 
absolute velocity ft/sec 
specific volume T UH 
relative velocity ft/sec 


no. of blades --- 


absolute flow angle degrees 
relative flow angle degrees 
angles defined by Figure 6-3 degrees 
angles defined by Figure 8-2 degrees 
"acute angle between y-axis and F, degrees 
acute angle between F, and m degrees 
stress psi 
rotational speed rad/sec 
increment === 
density lb/ft? 
degree of reaction cuc 
efficiency A 
Definition 


centrifugal or center 
center of mass 

gas force 

inner 

axial 


minimum 
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Subscripts Definition 


max maximum 

0 outer 

R rotor 

S stator, isentropic, stage 
ES y coordinates 


overscore vector 


General 

Of the quantities that were assumed in the analysis of the turbo- 
machines the most critical were the allowable peripheral speeds. These 
had an effect on the machine efficiencies, radial dimensions and the 
number of stages required. The choice of the allowable speeds was made 
for both the compressors and turbines on the premise that the blade and 
rotor stresses would not exceed safe limits. It is the objective in 
this chapter to investigate the parts of the machines where the stresses 
will be the greatest; namely, the first and last stages of Compressor 


No. l and the first stage of Turbine No. 1. 


Compressor Blades 





Figure 6-1 Forces Acting on a Compressor Rotor Blade 
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constant 


V = constant 
m 


From the first and last stage data given in Tables 4-2 and 4-5, 


the equations above, and Figure 4-3, the velocity diagrams for the 
r, center and outer radii have been determined as indicated in 


les 8-1 and 8-2. 








Inner Center Outer Units 


19.435 20525 21.065 in. 


616 604 593 








TABLE 8-2 


VELOCITY DIAGRAM DATA FOR THE LAST STAGE OF COMPRESSOR NO. 1 





Inner Center Outer Units 
E 19.435 20.10 20.765 in. 
SA 432 432 432 ft/sec 
V, 609 599 591 p 
E 0 0 0 a 
V $ 
u) 429 415 402 
N, 1262 1300 . 0333 " 
W, 869 920 964 ds 
W 
uy 1185 1225 1267 d 
W " 
u, 756 810 865 
U 1185 1275 1267 a 
t 
6 70 70.6 71.1 degrees 
G 2 60.2 62.0 63.4 2 
a | 0 0 0 , 
L, 44.8 43.8 43.0 " 
c 0.32 0782 0.82 in. 
Sp 0.79 0.82 0.85 a 
S 0.38 0.39 0.40 de 


Blade Profiles 
The profiles selected for the rotor and stator blades are NACA 


Airfoil 6309 and modified Joukowsky Airfoil 603, respectively. These 
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are illustrated in Figures 8-2 and 8-3 showing important dimensions. 
Table 8-3 lists important section properties. Angles Y and § 
depend on space/chord ratio, the velocity diagrams, and the type of 
profile. Actual values of Ö were obtained from de EL whereas 

Y | (15) | 
values for were obtained from Vavra . These are tabulated in 
Table 8-4. 

TABLE 8-3 


AIRFOIL SECTION PROPERTIES FOR A CHORD OF 0.82 INCHES 


NACA 6309 Joukousky 603 
Wen 0.042 0.0532 ae 
I. Ei 1.46 x 10? DD RENE S 
min AX 
=I O > o? 
max y 
TABLE 8-4 


VALUES OF Y AND § (FIGURES 8-2 AND 8-3) 


First Stage Last Stage 
inner Center Outer inner Center Outer 
P2 27.2 52 23.5 27.3 2503 23.8 
۱ 5 25.5 23.8 27.6 25.6 | 24.1 
JU 26.0 25.0 24.8 ---- — — 
x" 29.1 28.1 27.9 EN. An X 
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NACA 6309 





Location 


Figure 8-2 


x (in 


0.0928 
0.347 


0.473 


First Stage Rotor Blade 


140 





U 


3l 


e 


0.0492 
-0.041 


-0.034 


for Compressor No. 1 





Modified 
Joukousky Profile 603 





J 
۱ f D 
e > 
ES $ 
pn. 
= \ 
Coordinates of A: 
X 

x = 0,50l in. 
y = -0.1097 in. 
Figure 8-3 First Stage Stator Blade for Compressor No. 1 
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Centrifugal Blade Forces 


The maximum blade stress will occur at the blade root. A portion 
of this stress for the rotor blade will be due to the centrifugal force 


effects and can be calculated from the following equations: 


RR a en 
Ge pastas te (8-2) 


For constant chord blades of low alloy steel (P= 0.282 ns 


the centrifugal forces and stresses are as listed in Table 5-5. 


TABLE 8-5 
CENTRIFUGAL FORCES AND CORRESPONDING STRESSES FOR 


THE ROTOR BLADES OF COMPRESSOR NO, 1 


First Stage Last Stage Units 
133 733 rad/sec 
"eg 20.25 20.10 in. 
m 0.0193 0.0158 lb 
F 544 442 lb 
C T 
ue 12,950 10,500 psi 


Gas Bending Force 


The force acting on the blades due to the pressure differential 


and the change in momentum of the flow can be calculated by means of 
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the momentum theorem. 








(a) Rotor 


(b) Stator 
Figure 8-4 


Momentum Theorem Control Surfaces for the Blades 
of Compressor No. 1 


Referring to Figure 8-4 and using the dashed line as a control 


surface, the momentum equations become: 


For Rotor: 


my, سب‎ L. سس‎ mw 
WM + HEM = (RA - RAJA - o 
where a = 2 sh 
For Stator: 
DONE 7 tn E -3a 
- Hê V, + ŻY, = (RA, "RAM -E om 
where A 





Solving these equations for m requires knowledge of the number of 
blades, 7, and the pressure rise across the blades in addition to 
quantities previously determined. 


Za ITI blades (8-4) 


The pressure rise for the first stage is easily calculated from 


equation (1-6a) since T, and P, are known from Table 2-2 and Ah, from 


1 1 
Table 4-4. Although the degree of reaction, X, is defined in terms of 





enthalpies, it can be expressed in terms of pressure rise with little 


error. 


p hM STER, (415) 
ha-h, EE 


For the last stage, the quantities Ah, and T must be calculated. 


This can be done by equations (4-24) and (1-5a) since h,, U, T and Hs 


32 
are known. The pressure rise can then be calculated from equations 
(1-6a) and (4-15). 


Table 8-6 indicates values of Z, AP and Po: 
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TABLE 8-6 


GAS FORCES ACTING ON THE BLADES OF COMPRESSOR NO. 1 


First Stage Last Stage Units 
Rotor Stator Rotor Stator 
Z 155 326 144 324 Blades 
EP 2572 4.9 29.8 6.1 psi 
E 36.0 7.92 35.9 cue lbs 
8 26,2 66.7 25.2 --- degrees 


The gas bending stresses have been calculated for three locations 
on the rotor blades and one on the first stage stator blade. 


Referring to Figures 8-2 and 8-3, the stresses can be expressed as: 





ES ۱۳ . ۹ (8-5) 
E JE 
x 4 
where W — in: E, sin $ 
B p 
9 pe = Fq Cos Y 


eO بت‎ 2i for rotor 


b 
۱ 


= (e — 90) + Y! for stator 
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TABLE 5-7 


BLADE ROOT STRESSES DUE TO GAS FORCES 


First Stage Last Stage First Stage Units 
م‎ S -2.4 +4 .6 degrees 
H. 22085 298 -6.42 in - Ib 
SØ (0.62 1.0 +0.52 in - lb 
Es da -99 , 000 -60, 300 +12,850 psi 
( Tod +82 ,350 +66 , 600 == psi 
( Fade +68 , 0 +55,800 -..- psi 


Applying the centrifugal force stresses to the gas force stresses, 
a maximum stress of 95,300 psi occurs at point B of the first stage 
rotor blade. This stress is higher than normal long-life design criter- 
ia and would preferably be reduced to 50,000 psi or less. Lower 
stresses can be obtained by modifying the blades. Instead of using a 
9% thick blade at all sections, a tapered blade 12% thick at the root 
would reduce the stress to 53,700 psi. If, in addition, the chord is 
increased by 10%, the stress is further reduced to about 40,400 psi. 

There are many low-alloy, high-strength steels that would be suit- 
able for the compressor blades such as AISI Nos. 4140, 4340, 6150, 
8640 and 8740. 

As one might expect, the stresses in the stator are small compared 
to the rotor stresses since most of the pressure rise is in the rotor 
blades. Providing there is no serious vibration problem, no major 


design changes would have to be made to permit a mean peripheral speed 
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of 1235 ft/sec. 





Turbine Blades 


en 


I 


ty 


DT 


| 


۷ 


Figure 8-5 Forces Acting on a Turbine Rotor Blade 


For free-vortex flow and equal work output at all radii: 
rv = constant 


= constant 
m 


From the first stage data given in Table 5-1, the equations above, 
and Figure 5-3, the velocity diagrams for the inner, center and outer 


radii have been determined as indicated in Table 8-5. 
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Also: 


N 


7000 rpm 
ERC ODIN, 
85 


71 
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Units 


in. 


ft/sec 








Rotor Blades 

The root, center and tip profiles of the first-stage rotor blade 
are illustrated in Figure 8-6. This is a tapered blade in which the 
centroids of each section lie on a radial line. The maximum blade 
stresses occur at the blade root, hence only these section properties 
have been determined. As in the compressor, the stress is due to a 
combination of gas bending force and centrifugal force effects. 

Centrifugal Force 

The volume and center of mass of the tapered blade having root, 
center and tip sections as shown in Figure 8-6 can be determined from 


the following plot: 


0 
® 


2) 


Section Area (in 


0 2.3 526 
Distance from Blade Root (in.) 


Blade Volne de= 0.21 + 0.48 meer 1.82 in^ 
T - 10.94 + 2.30 = 13.24 in 


CG 


For alloy 5-816: ( (E = 0.310 In) 


= 
u 


0.564 lb per blade 


F = 10,500 1b per blade 


S 


21,900 psi 
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Tip 


Center 


Root 


Figure 8-6 


=] 


150 


Arca = 0 TIT 


Spacing = 1,19 in. 


Area = 0,345 in? 


Spacing =» 1.0 in. 


Area = 0.48 in^ 


Spacing s 0.81 in. 


IST . € 0.014 0858 
min 


y “max 


>x 


First Stage Rotor Blade Profiles for Turbine No. 1 





Gas Force 





Before the gas force can be calculated by means of the momentum 
theorem, the pressure differential betwéen the inlet and exit to the 
rotor blade must be calculated. 


Referring to Figures 5-4 and 5-5, the following relations are 











apparent: 
2 2 
= L = ۱ ۰ ۸ ۱ 
(Ah), = L2 p* 29 Pao (8-6) 
ha = h, + (1-%) 5 (5-12) 
2 a 
h, = h, - e | 
(ah) = (ah - بط‎ é on 
Also: 
| QVE 
( N), Th pli - E ou 
۳ Der Dur ۲ ون‎ (1-5) 
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* 
AR eE Ane (5-10a) 
ls 
(ahy با‎ + 
A ZNA K K 2 5-10b 
5) ( S 293 ( ) 
Combining equations (8-8), (5-10a) and (5-10b): 
2 2 
93 Ah MS Ve 
2 < کم‎ - > )8-9( 
( s) s ETE p^ 29) 
Also: 
Y 
9 PA € 
(ahs) Lh Cp e p (1-7a) 
E 
From these equations, Tables 2-2, 8-92 and Chapter V, the pressure 
rises are: 
Fi z Po LED psi 


Using the dashed line of Figure 8-7 as a control surface, the 
momentum theorem becomes: 


Ee 


ms 








m. vs __ DIS 
PES No + -F W, = (RA, - RA) 7- Ë 
where Ay = A, = sh 
E = 60.5 7 + 7200 
MEE ها‎ © — 499° 


152 





50216 20304 
SUTqiN] 9 103] 5 
3[1J014 2]2Z0N a3urqunj 9-9 31N314 ] 013109 Wa3109U] ۲۲3 0 L-8 21۳327 آر‎ 


7۰] = Sutoeds 


I 


e 


ix 
É 
vo] 








153 








The fibers most heavily stressed in the root cross section are 


indicated by A and B in Figure 8-6. Axes x and y are principal axes. 





pee Mox (8-5) 
I T, 


(16) 


A total stress of 29,000 psi exists at A. Jennings states that 
the blade root temperature will normally be 175-300ºF lower than the 
gas stagnation temperature or a temperature of about 1100ºF here. The 
stress-to-rupture in 100,000 hours for S-816 alloy is 47,000 psi at 
1100ºF. Therefore, the rotor blade stresses are not excessive for the 
blading of Figure 8-6 providing no serious vibration problem exists. 
stator Blades 
There is no attempt here to do a detailed stress analysis on the 
nozzles. Figure 8-8 illustrates a typical profile for the root section. 


The gas force has been calculated to be 200 lbs at © = 25.6°. The 


minimum section modulus, "E , of the profile is approximately 17.6 x 


max 
107? I and can be increased by increasing the % thickness and/or 


chord lenzth. 
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( = Mx Ymax = 30,000 psi 
max T 
x 


Since the stator blade temperatures will be higher than the rotor 
blades, a maximum stress of less than 25,000 psi is desirable and this 


can easily be achieved by slight modification of the blade. 


Turbine Disks 


o) 


(b) 





Figure 6-9 Turbine Disk 
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The type of disk being considered is illustrated in Figure 8-9. 
All dimensions except the disk width, b, are known from blade data and 


practical considerations. The value of b depends on stress conditions 





due to the centrifugal forces. The stresses in the rim (r,<r<r,) are 
substantially smaller than the blade root stresses so they need not be 
calculated here. Since we are interested in the stresses and dimensions 
of the remainder of the disk, the blade and disk assembly can be replac- 
ed by the hub portion only (O<r<r,) with an external load due to 
centrifugal forces of the mass of the blades and the rim (Figure 8-9(b)). 
Because of the irregular shape of the rim it is divided into sec- 


tions (Figure 8-10) for the purpose of determining centrifugal forces. 


p 


Figure 8-10 Approximate Cross Section of Turbine Disk Rim 
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TABLE 8-9 


CENTRIFUGAL LOADING ON THE TURBINE DISKS 


Item Material P (b/in”) “CG (in) m(lb/in at A) q(1b/in at A) 
1 Fir Tree: S-816 
Rim: 16-25-6 0.300 10.145 0.957 13,600 
2 16-25-6 0.291 9.170 0.150 1,920 
3 16-25-6 0.291 8.845 0.266 3,270 
4 16-25-6 0.291 8.995 0.126 1,580 


q,, = 20,370 
m 


Loading at A due to blades; 


blade: joe cic) 16,350 ib/in at A 
2x A 
Total loading: 
q(total) = 16,350 + 20,370 = 36,720 lb/in at A 


For a constant stress disk: 





pd (mir 
b= be *9" (8-10) 


Since the disk temperature will be approximately 1000ºF at A and 
500ºF at the center, an allowable stress of 25,000 psi is assumed, 
based on 100,000 hour stress-to-rupture data of 16-25-6 alloy. 


ai = Sexo — ¡41 


E = in. 
q ES EO 
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WIDTHS OF CONSTANT STRESS DISK 


b (in.) 


3.15 


3.12 


3.02 


2.87 





Figure 8-11 illustrates the first and second stage rotor and 
nozzle assemblies. 

It appears from the foregoing stress calculations that the peri- 
pheral speed that was assumed is not excessive. There is, of course, 
the possibility of much higher blade stresses if the unsteady gas force 
frequency is close to the blade natural frequencies. This problem will 
not be investigated here since normally these effects can be minimized 
by modification of the blading rather than major redesign of the 


turbine. 
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CHAPTER IX 


SUMMARY AND CONCLUSIONS 


In this report, a limited design study was performed for each of 
the major components of a 50,000 HP closed-cycle gas turbine power 
plant employing helium a$ the working fluid and a nuclear reactor as 
the heat source. The descriptive data tabulated in Table 9-1 should be 
considered representative rather than optimum for a cycle having a maxi- 
mum pressure of 1000 psia and a maximum temperature of 1350°F. It is 
believed that comparable sizes and performance would be obtained if 
different choices are made for such variables as compressor and turbine 
rpm, type of blading, ratio of blade spacing to chord, type of heat 
exchanger surfaces, and reactor cooling channel arrangement, 

The thermal efficiency of the cycle is in the order of 33-36% 
depending on the accuracy of the machine efficiencies and pressure 
drops. For temperatures of 1450 - 1500°F it is estimated that the 
cycle efficiency would increase to about 40%. The size of the compon- 
ents would be reduced substantially from those herein if higher pres- 
sures were employed, 

Although this design study was based on a ship propulsion plant, 
no mention has been made as to its operational requirements. One of 
these requirements is astern power. It may be possible to use an astern 
gas turbine in an arrangement similar to that of a conventional steam 
propulsion plant. 

Another operational problem of critical importance concerns the 


control of the cycle. The primary control consists of raising and 
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^ FOR THE MAJOR COMPONENTS 


| Max. Blade Tip Max. Blade 
Unit RPM Diameter (in.) Height (in.) Efficiency 


Compressor #1 7000 1 1.6 0.867 
Compressor #2 7000 40.9 no 0.840 
۱ à 


Turbine #1 7000 | 6.2 0.858 


Turbine #2 5500 l 7.1 0.834 





Overall Outer Tube 








Length Diameter Length ID. No. of AI B 
(f) (ft) (ft) in.) Tubes i 
Regenerator 45 1525 2557 0.30 10,330 0.0194 
Precooler 18 4.6 10.2 0.40 2,690 0.0158% 
Intercooler 18 30 9.9 0.40 1, 700 0.0062 
Reactor Core 1.565 ED 1205 1.08 2,540 0.0054 


* Includes the low pressure side of the regenerator. 
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lowering the upper pressure level by means of an auxiliary system con- 
taining high pressure and low pressure helium accumulators. The speed 
with which this control functions will have an effect on ship maneuver- 
ability. Due to the large inventory of helium in the cycle this method 
of control may be sluggish compared to that of conventional steam 
plants. One important characteristic of this type of cycle is its good 
efficiency over a wide range of power. 

Since helium is unavailable at sea except for a limited reserve 
supply, there must be practically zero leakage. This is very difficult 
to achieve when rotating machinery is involved. The main leakage prob- 
lem occurs at Turbine No. 2 where the shaft leaves the casing. The 
remainder of the cycle can be enclosed in a leakproof casing. Bearing 
lubrication then becomes a problem. If oil is used as the lubricant, a 
rather complex system of separating helium from the oil is required in 
order that the helium can be returned for further use. It is possible 
to design the bearings so that high pressure helium is the lubricant 
while the machines are operating. With this arrangement, an auxiliary 
oil system would be required for starting. 

Long half-life fission products must be prevented from being car- 
ried to all parts of the cycle by proper design of the cooling channels 
and filtering system. Then, only the reactor would have to be shield- 
ed. 

These are only some of the problems that would have to be thor- 
oughly investigated before a meaningful appraisal of the closed-cycle 
helium gas turbine plant can be made, However, on the basis of the 
limited investigation of this report, such a cycle has definite advan- 
tages and is worthy of further study. 
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